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This  paper  presents  the  results  of  a study  to  develop  a method  for 
predicting  various  stochastic  characteristics  of  wave  groups  in  random 
seas  where  the  water  depth  is  sufficiently  deep.  Wave  groups  are 
defined  in  the  present  study  as  a sequence  of  at  least  two  high  waves 
having  nearly  equal  periods.  Formulae  to  statistically  predict,  for  a 
given  wave  spectrum,  (i)  time  duration  of  a wave  group,  (ii)  number  of 
waves  in  a group,  and  (iii)  frequency  of  occurrence  of  wave  groups  in  a 
specified  time  are  newly  derived.  The  probability  density  functions  of 
time  duration  of  wave  group  as  well  as  time  interval  between  successive 
wave  groups  are  analytically  developed  based  on  the  concept  of  a narrow- 
banded  Gaussian  random  process.  Renewal  theory  is  applied  to  the 
probability  density  function  of  time  interval  between  wave  groups  for 
evaluating  the  frequency  of  occurrence  of  wave  groups  in  a specified 
time.  The  results  of  numerical  computations  carried  out  using  wave  data 
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measured  in  the  North  Sea  show  that  the  average  time  duration  of  a wave 
group  computed  by  the  present  theory  agrees  reasonably  well  with  the 
observed  data.  A.lso,  the  frequency  of  occurrence  of  wave  groups 
predicted  by  the  present  theory  agrees  well  with  those  observed,  while 
the  frequency  of  occurrence  computed  by  applying  the  formula  currently 
used  is  substantially  greater  than  observed  data. 
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CHAPTER  I 
INTRODUCTION 


A unique  phenomenon  associated  with  the  wind  generated  waves  in  the 
ocean  is  a sequence  of  high  waves  having  nearly  equal  periods  commonly 
known  as  wave  groups.  Two  examples  of  wave  groups  observed  in  the  ocean 
are  shown  in  Figure  1.1.  Figure  1.1(a)  is  taken  from  Rye  (1974)  in 
which  he  shows  group  waves  recorded  in  the  North  Sea  during  a storm  of 
® f i cant  wave  height  of  10.0  meters,  while  Figure  1.1(b)  shows  wave 
groups  recorded  in  severe  seas  by  weather  ship  at  Station  K in  North 
Atlantic. 

It  has  been  known  that  wave  groups  often  cause  problems  for  safety 
of  marine  systems  when  the  period  of  the  individual  waves  in  the  group 
is  closed  to  the  marine  system's  natural  motion  period.  This  is  not 
because  wave  heights  are  exceptionally  large  but  occurs  primarily 
because  of  motion  augmentation  due  to  the  resonance  with  the  waves 
which,  in  turn,  can  induce  capsizing.  As  another  example,  a moored 
marine  system  tends  to  respond  to  successive  high  waves  which  induce  a 
slow  drift  oscillation  of  the  system  resulting  in  large  forces  on  the 
mooring  lines. 

Since  the  wave  group  phenomenon  is  extremely  interesting  and  impor- 
tant for  the  consideration  of  marine  systems,  many  studies  have  been 
carr^e<^  ou*-  on  the  various  phases  of  the  phenomenon.  One  area  of  con- 
siderable interest  is  the  stochastic  analysis  of  wave  groups  in  random 
seas . 
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Examples  of  wave  groups  observed  at  sea 
(Example  (a)  from  Rye  (1974)). 


Figure  1.1 
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Several  studies  have  been  carried  out,  to  date,  on  the  stochastic 
analysis  of  wave  groups.  However,  in  almost  all  these  studies,  the 
length  of  time  that  the  wave  group  persists  (Ta+  in  Figure  1.2),  the 
number  of  waves  involved  in  a group,  and  the  time  interval  between 
successive  occurrences  of  wave  groups  (ta  in  Figure  1.2)  are  all  average 
values  and  no  precise  probabilistic  information  on  these  individual 
quantities  has  been  provided. 

Furthermore,  extreme  care  must  be  taken  in  applying  the  currently 
available  methods  for  evaluating  the  mean  number  of  waves  and  mean  time 
of  duration  of  wave  groups  due  to  the  following  reason.  It  has  been 
customary,  to  date,  to  consider  the  exceedance  of  the  envelope  above  a 
certain  level  to  identify  a wave  group.  However,  this  is  not  correct; 
if  the  time  duration  shown  in  Figure  1.2  is  relatively  short,  there 
may  be  one  wave  crest  (or  no  wave  crest)  in  time  xa+,  which  obviously 
does  not  constitute  a wave  group  even  though  the  envelope  exceeds  the 
specified  level.  In  order  to  elaborate  this  statement.  Figure  1.3  shows 
a pictorial  sketch  of  the  time  interval  between  successive  upcrossings 
of  the  envelope  of  a random  process.  Suppose  the  upcrossing  of  the 
envelope  at  point  B in  the  Figure  1.3  is  not  associated  with  wave 
groups,  then  the  time  interval  AC  (instead  of  AB  and  BC)  should  be 
considered  as  time  interval  between  wave  groups. 

It  is  the  purpose  of  this  study  to  develop  an  analytical  method  for 
evaluating  various  stochastic  characteristics  of  wave  groups  in  random 
seas  from  a given  wave  spectrum.  The  water  depth  is  considered  to  be 
sufficiently  deep  in  the  present  study.  The  analytical  method  includes 


the  evaluation  of 
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Figure  1.2  Level  crossing  of  the  envelope  of  a random  process 
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Figure  1.3  Definition  of  time  interval  between 


successive  wave  groups,! 
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(a)  The  probability  density  function  of  the  time  duration  of 
wave  group, 

(b)  The  probability  of  occurrence  of  a specified  number  of  waves 
in  a group, 

(c)  The  probability  density  function  of  the  time  interval 
between  successive  wave  groups,  and 

(d)  Frequency  of  occurrence  of  wave  groups  in  a given  sea. 

In  the  derivation  of  the  probability  density  functions  of  this 
study,  the  wave  group  phenomenon  is  treated  as  a level  crossing  problem 
associated  with  the  envelope  of  a narrow  band  Gaussian  random  process 
(R(t)  in  Figure  1.2),  and  then  the  probability  density  functions  are 
modified  taking  the  condition  required  for  the  wave  groups  into  con- 
sideration. Wave  groups  are  defined  in  the  present  study  as  a sequence 
of  at  least  two  high  waves  exceeding  a specified  level.  For  evaluating 
the  frequency  of  occurrence  of  wave  groups  in  a specified  time,  renewal 
theory  is  applied  to  the  probability  density  function  of  the  time  inter- 
val between  wave  groups. 

The  average  time  duration  associated  with  wave  groups  computed  by 
the  newly  developed  theory  and  the  predicted  frequency  of  occurrence  of 
wave  groups  agree  well  with  the  data  observed  in  the  North  Sea. 

This  paper  consists  of  three  chapters  followed  by  conclusions  and 
eleven  appendices.  Chapter  II  reviews  the  currently  available  litera- 
ture on  the  stochastic  characteristics  of  wave  groups  in  random  seas 
where  the  water  depth  is  sufficiently  deep. 

In  Chapter  III,  the  time  duration  of  wave  group  and  the  number  of 
waves  in  a wave  group  are  considered.  The  probability  density  function 
applicable  for  the  length  of  time  that  the  wave  group  persists  is 
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analytically  developed.  Formula  to  evaluate  the  probability  of  the 
number  of  waves  in  a group  is  derived  based  on  the  probability  density 
function  of  time  duration.  The  average  time  duration  associated  with 
wave  groups  computed  by  the  newly  developed  theory  was  compared  with  the 
data  observed  in  the  North  Sea. 

In  Chapter  IV,  a method  to  estimate  the  frequency  of  occurrence  of 
wave  groups  is  discussed.  The  probability  density  function  of  time 
intervals  between  successive  wave  groups  is  derived  and  then  a method  to 
evaluate  the  frequency  of  occurrence  of  wave  groups  in  a specified  time 
is  developed  by  applying  renewal  theory.  Measured  data  in  the  North  Sea 
and  the  predicted  frequency  of  occurrence  of  wave  groups  are  compared. 

In  Chapter  V,  the  findings  and  conclusions  obtained  from  the 
results  of  the  present  study  are  summarized. 

In  the  Appendices,  additional  derivations  to  clarify  or  to  prove 
formulae  used  in  the  main  text  are  presented. 


CHAPTER  II 

REVIEW  OF  THE  LITERATURE 


This  chapter  reviews  the  currently  available  literature  on  the 
stochastic  characteristics  of  wave  groups  in  random  seas. 

Two  different  approaches  have  been  applied  to  the  stochastic 
analysis  of  the  wave  group  phenomenon.  One  approach  is  to  consider  wave 
groups  as  a level  crossing  problem  associated  with  the  envelope  of  random 
processes.  The  other  approach  regards  the  sequence  of  wave  heights  as  a 
Markov  chain  with  a correlation  only  between  two  successive  waves. 

The  first  approach  was  used  by  Rice  (1945,  1958)  in  his  pioneering 
work  on  the  mathematical  analysis  of  a Gaussian  noise  and  was  first 
applied  to  the  stochastic  analysis  of  wave  groups  in  random  seas  by 
Longuet-Higgins  (1957). 

Longuet-Higgins  (1957,  1984)  derived  formulae  to  evaluate  the  mean 
length  of  the  run  (i.e.,  the  mean  number  of  waves  included  in  a wave 
group),  the  mean  time  interval  between  successive  upcrossings  of  a 
specified  level  by  an  envelope,  and  the  mean  time  duration  of  a wave 
group,  in  terms  of  spectral  moments  and  the  band  width  parameter  of  the 
frequency  spectrum.  Assuming  that  the  upcrossings  by  the  wave  envelope 
of  a certain  level  were  related  with  wave  groups  and  followed  a Poisson 
random  process,  the  intervals  between  successive  wave  groups  were  shown 
to  be  exponentially  distributed.  Furthermore,  he  stated  that  the  two 
approaches  (wave  envelope  approach  and  the  Markov  chain  approach)  led  to 
almost  equivalent  results  and  both  could  be  applied  only  to  sufficiently 
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narrow  band  random  processes,  and  to  data  in  which  higher  and  lower 
frequencies  were  filtered  out. 

The  mean  number  of  waves  included  in  a run  of  high  waves  and  the 
average  number  of  waves  between  successive  occurrences  of  wave  groups 
were  given  by  asymptotic  formulae  derived  by  Ewing  (1973).  The  frequency 
spectrum  was  assumed  to  be  symmetric  about  its  mean  frequency  and  narrow 
banded.  Also  presented  were  the  asymptotic  formulae  as  functions  of  the 
spectral  peakedness  parameter,  which  represented  the  degree  of  concentra- 
tion of  the  frequency  spectrum  about  its  mean  frequency.  When  the  spec- 
tral peakedness  parameter  had  large  values  (very  narrow  banded  frequency 
spectra),  the  predicted  mean  values  by  Ewing's  asymptotic  formulae  agreed 
well  with  data  from  numerical  wave  simulation  experiments.  However, 
Chakrabarti  et  al.  (1974)  commented  that  the  mean  length  of  the  run  of 
high  waves  between  successive  occurrences  of  wave  groups  was  smaller  than 
the  observed  values  in  wave  data  from  a North  Atlantic  storm. 

Nolte  and  Hsu  (197  2)  studied  the  average  time  duration  of  the  wave 
envelope  exceedances  over  a specified  level  (denoted  by  Ta+  in 
Figure  1.2),  the  average  frequency  of  occurrences  of  the  upcrossings  by  a 
wave  envelope  of  a high  level,  and  the  average  time  between  successive 
level  upcrossings  by  an  envelope  (^  in  Figure  1.2).  The  wave  envelope 
upcrossings  of  a certain  level  were  considered  to  be  related  with  wave 
groups.  By  assuming  that  the  wave  group  occurrences  were  a Poisson 
random  process,  a formula  was  derived  to  evaluate  the  probability  of  a 
specified  number  of  wave  groups  in  a certain  time  interval. 

Goda  (1976)  assumed  that  the  successive  wave  heights  were  statis- 
tically independent.  By  applying  the  theory  of  runs  (Fisz,  1963),  he 
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calculated  the  probability  that  a sequence  of  wave  heights  higher  than  a 
certain  level  (i.e.,  a wave  group)  contained  a specified  number  of  wave 
heights.  Using  the  envelope  theory,  the  mean  time  duration  of  the 
envelope  persistence  above  a specified  level  was  derived.  Also  developed 
was  a formula  for  evaluating  the  mean  length  of  the  run  of  high  waves  as 
a function  of  the  spectral  peakedness  parameter.  The  predicted  mean 
number  of  waves  in  a wave  group  was  compared  with  that  which  was  observed 
in  field  wave  data  and  numerical  wave  simulation  data,  and  he  concluded 
that  the  degree  of  the  spectral  peakedness  influenced  the  mean  number  of 
waves  in  a wave  group.  However,  Wilson  and  Baird  (197  2)  found  that  the 
mean  value  of  the  number  of  waves  included  in  a wave  group  evaluated  by 
using  Goda  s formula  was  smaller  than  the  observed  value  in  wave  data 
obtained  in  the  ocean  off  Nova  Scotia. 

Rye  (1979,  1981)  found  that  the  degree  of  the  spectral  peakedness 
was  a significant  factor  for  the  wave  group  formation  and  noted  that  the 
correlation  coefficient  of  two  consecutive  wave  heights  was  strongly 
related  with  the  spectral  peakedness  parameter.  Furthermore,  Rye  (1974) 
observed  in  wave  data  from  the  North  Sea,  that  the  wave  group  formation 
appeared  to  be  more  pronounced  during  the  growing  stage  of  wind  generated 
seas  than  during  the  decaying  stage. 

Mase  et  al.  (1983)  calculated  the  mean  length  of  wave  groups  (i.e., 
the  mean  number  of  waves  in  a wave  group)  and  the  mean  number  of  waves 
between  consecutive  wave  groups  by  using  numerical  wave  simulation  and 
field  data.  Both  average  values  from  the  two  types  of  data  were  compared 
in  order  to  study  the  degree  of  wave  grouping. 
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Elgar  et  al.  (1984)  argued  that  the  mean  length  of  wave  groups  could 
be  generally  related  with  the  spectral  peakedness  parameter.  They 
calculated  values  of  the  mean  length  of  wave  groups  from  numerical  wave 
simulation  data,  plotted  them  against  the  spectral  peakedness  parameter, 
and  found  that  the  spectral  peakedness  parameter  and  the  mean  length  of 
the  wave  groups  are  well  correlated  for  broad  banded  frequency  spectra. 

The  second  approach  that  regards  a sequence  of  wave  heights  as  a 
Markov  chain  was  used  first  by  Kimura  (1980).  The  correlation 
coefficient  of  consecutive  waves  was  evaluated  from  numerical  wave 
simulation  data,  and  thereby  estimates  of  the  mean  number  of  waves  in  a 
wave  group  and  the  mean  number  of  waves  between  successive  wave  groups 
were  obtained. 

Battjes  and  Vledder  (1984)  estimated  the  correlation  coefficient  of 
two  consecutive  wave  heights  by  using  a spectral  shape  parameter  based  on 
the  theory  of  Gaussian  processes.  Their  results  for  the  mean  length  of 
the  run  and  for  the  mean  number  of  waves  between  successive  wave  groups 
agreed  with  those  derived  by  Kimura. 

Arhan  and  Ezraty  (1978)  developed  a formula  for  the  joint 
probability  density  function  of  two  consecutive  wave  heights  that  can  be 
used  to  calculate,  by  the  means  of  Markov  chain  theory,  the  average 
number  of  waves  in  a wave  group. 

As  this  brief  review  indicates,  several  studies  have  been  conducted 
for  the  stochastic  characteristics  of  wave  groups  in  random  seas. 
However,  in  almost  all  these  studies,  the  number  of  waves  included  in  a 
wave  group  and  the  time  duration  of  a wave  group  over  a certain  level  are 
average  values,  and  no  precise  probabilistic  formulae  for  these 
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individual  quantities  have  been  provided.  The  analysis  of  wave  records 
has  indicated  that  the  currently  available  theories  underpredicted  the 
mean  length  and  the  mean  time  duration  of  wave  groups.  Furthermore,  most 
of  the  studies  on  wave  groups  assumed  that  the  envelope  upcrossings  of  a 
specific  level  were  related  with  wave  groups  and  followed  a Poisson 
random  process.  However,  it  is  highly  desirable  to  verify  the  validity 
of  this  assumption  for  evaluating  the  frequency  of  occurrence  of  wave 
groups  in  random  seas. 


CHAPTER  III 

TIME  DURATION  OF  AND  NUMBER  OF  WAVES  IN  A WAVE  GROUP 


Probability  Distribution  of  Time  Duration  of  Envelope 
Exceeding  a Certain  Level 

As  stated  in  the  Introduction,  an  exceedance  of  the  envelope  above 
a specified  level  a does  not  necessarily  mean  a wave  group  unless  two  or 
more  wave  crests  are  compromised  in  the  exceedance.  Therefore,  we  first 
derive  the  probability  density  function  of  the  time  duration,  xa+ 
associated  with  the  envelope  exceeding  a level  a;  then  this  probability 
density  function  will  be  truncated  taking  into  consideration  the 
condition  required  for  the  existence  of  a wave  group. 

In  the  development  of  the  probability  density  function  of  time 
duration  of  the  envelope  exceeding  a specified  level  a,  waves  are 
assumed  to  be  a Gaussian  random  process  with  a narrow-band  sprectrum. 
Then,  following  Rice's  approach  (1945),  the  probability  that  the 
envelope  of  a random  process  x(t)  exceeds  a specified  level  a at  time 
with  velocity  during  an  upward  crossing  and  then  crosses  that  same 
level  downward  at  time  t2  with  velocity  R2  is  given  approximately  by  the 
following  formula: 

i" Upward  crossing  of  a level  a at  time  t, 

I 1 

Pr  ■ with  velocity  R^  followed  by  a downward 
_ crossing  at  time  t2  with  velocity  R2 


]=  v 

i v 
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'i  t t ii 

R1R2  f(o,R1>a,R2)  dR1  dR2 


Rj  f(a,R  ) dRL 


o < R^  < °° 

-“>  < R2  < o (3.1) 


Here,  the  numerator,  ^a+(Ta+)>  represents  the  expected  number  (per 

unit  time)  of  envelope  upward  crossings  with  velocity  R followed  by  the 

downward  crossings  with  velocity  R_,  while  the  denominator,  N 

L a+* 

represents  the  expected  number  (per  unit  time)  of  envelope  upward 
crossings  of  the  level  a with  velocity  R^ . The  function  f(a,  Rp  a,  R,) 
is  the  joint  probability  density  function  of  displacement  R and 

t 

velocity  R at  times  t^  and  t2  with  R^  = R2  = 01  in  this  case,  while 

1 

f(“,  is  the  joint  probability  density  of  displacement  R and 

1 

velocity  R of  the  envelope  at  time  t ^ with  Rj  = a. 

The  probability  given  in  Eq.  (3.1)  may  be  interpreted  as  the  time 
that  the  envelope  spends  above  the  level  a.  In  other  words,  by  letting 

tl  ~ t2  Ta+*  (3.1)  is  equivalent  to  the  probability  density 

function  of  for  a specified  level  a. 

Integration  of  the  joint  probability  density  function  involved  in 
the  numerator  of  Eq.  (3.1)  cannot  be  carried  out  directly.  It  is 
necessary  to  first  obtain  the  joint  probability  density  function  of  the 
wave  profile  x(t)  and  velocity  x(t)  at  times  tl  and  t2,  and  then  to 
transform  this  joint  probability  density  function  of  displacement  and 
velocity  to  that  of  the  wave  envelope  and  its  velocity. 

Assuming  a narrow-band  Gaussian  random  process,  the  wave  profile 
x(t)  can  be  written  as 
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x(t) 


where 


“ ^ cn  cos(unt  + e^)  = xc(t)  cos  0)  t - x (t)  sin  u>  t (3.2) 


n=l 


xc(t) 


oo 


cos  {(u»  - w)t  - e } 
n n 


n=l 


(3.3) 


co 

,(t)  = / c sin  { (oj  - ZJ)t  - e } 

> t—i  n n n 


n=l 


“ = mean  frequency  of  wave  spectrum 
Let  the  random  variable  Xc^  represent  the  value  of  x£  (t)  at  time 
tl>  and  the  random  variable  Xgl  represent  the  value  of  x (t)  at  time  t, 
etc.  Since  we  take  sine  and  cosine  components  for  both  displacement  and 
velocity  at  times  tj  and  t2,  we  have  a set  of  eight  random  variables 

' t ! , 

cl’  Xcl’  Xc2*  Xc2’  Xsl’  Xsl’  Xs2’  Xs2^  which  composes  a random  vector 
X.  Here,  each  element  of  X obeys  the  normal  probability  distribution, 
and  we  may  write  the  joint  normal  probability  density  function  in  the 
following  form: 


f(X)  = l— 

(2  IT) 


= 1 


~rh 


-{x'S-'x 


(3.4) 


X’  = (x 


cl’  xcl’  XC2’  xc2’  Xsl’  xsl’  Xs2’ 


s2  ’ 


where 
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and  Z is  the  covariance  matrix  given  by 


Z = 


where 


U = 
o 


m 


o 


S(u))  dm 


U 


2 


v 

o 


V 


1 


V 


2 


n 

0 


n 


l 


oo 


00 


00 


00 


(m-m)  ^ 


S(m) 


du) 


m S ( (o)  dm  / 


S ( oj)  dm 


S(m)  cos(u>-u))t  da) 


S(m)  (m-m)  sin  (m  - m)x  dm 


S(m) 


_ 2 _ 

(m-m)  cos(m-m)  x dm 


S(m)  sin(m-m)  x dm 


(3.5) 


(3.6) 


S(m)  (m-m)  cos(m-m)  x dm 
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T 


t 


2 


t 


S(oj)  = wave  specrum 


The  detailed  derivation  of  the  elements  of  the  covariance  matrix  is 
given  in  Appendix  A. 

As  can  be  seen  in  Eqs.  (3.4)  through  (3.6),  the  joint  probability 
density  function  f(X)  can  be  evaluated  from  a given  wave  spectrum.  In 
practice,  however,  the  inverse  operation  of  the  covariance  matrix 
involved  in  Eq.  (3.4)  is  extremely  difficult  to  perform.  Even  though 
the  inverse  can  be  achieved,  the  follow  up  transformation  to  the  joint 
probability  density  function  of  the  envelope  and  its  velocity  is  not 


One  way  to  overcome  this  difficulty  is  to  decompose  a given 
spectrum  into  two  parts,  each  part  being  symmetric  about  its  mean 
frequency.  One  contains  primarily  the  lower  frequency  energy  components 
of  the  wave  spectrum,  while  the  other  contains  the  higher  frequencies  of 
the  spectrum  as  illustrated  in  Figure  3.1.  Each  part  of  the  spectrum 


feasible. 


can  be  expressed  in  the  form  of 


(3.7) 


SPECTRAL  DENSITY 
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Figure  3.1  Decomposition  of  wave  spectrum  into  two  symmetric- 
shape  spectra. 
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where  i = 1,  2 

“mi  = nodal  frequency  of  S.(oj) 
ai  - value  of  S^(u)  at  the  modal  frequency 
moi  = area  under  the  spectrum  S^rn). 

The  sum  of  the  two  areas  mQl  and  mQ2  is  equal  to  the  area  under  the 
original  spectrum  S(m)  which,  in  turn  is  equal  to  the  wave  variance, 

V Since  the  variance  of  the  velocity,  in  the  covariance  matrix 
also  plays  a significant  role  in  developing  the  probability  density 
function  of  time  duration  Ta+,  the  sum  of  the  variances  of  the  two 
decomposed  spectra  is  maintained  constant  and  is  equal  to  the  variance 

U2  of  the  original  spectrum.  Under  these  conditions,  the  parameters 

V “mi’  and  moi’  are  determined  numerically  through  a nonlinear  least 
squares  fitting  technique  such  that  the  difference  between  the  shape  of 
the  original  spectrum  and  the  sum  of  the  two  symmetric  spectra, 
expressed  by  | S (to)  - 2 S±(oj)  | 2,  is  minimal.  A detailed  discussion  on 
the  decomposition  of  the  spectrum  is  given  in  Appendix  B. 

Figure  3.2  shows  an  example  of  wave  spectrum  expressed  by  the  sum 
of  two  symmetric  spectra.  The  original  spectrum  represents  the  two- 
parameter  wave  spectral  formulation  for  a significant  wave  height  of 
10.0  m with  a modal  frequency  of  0.42  rps.  Although  there  is  some 
discrepancy  in  the  magnitude  of  spectral  densities  between  the  original 
and  the  combined  symmetric  spectra  at  frequencies  ranging  from  0.60  and 
0.85  rps,  the  values  of  moments  which  are  significant  for  the  analytical 

solution  of  wave  groups  are  very  close.  That  is,  y = 6.25  m2  and 

2 2 a 
U2  = 0.317  m /sec  of  the  combined  spectra  are  chosen  to  be  equal  to 

those  of  the  original  spectrum.  The  computed  moments  mj  and  m2  are  3.36 


SPECTRAL  DENSITY  IN  M2-SEC 
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Figure  3.2  Comparison  of  two-parameter  wave  spectra  ( significant  wave 
height  10.0  m. , modal  frequency  0.42  r.p.s.)  and  sum  of  two 
symmetric-shape  spectra. 
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2 o o 

m /sec  and  2.22  m /sec  , respectively,  for  the  combined  spectra  as 
compared  with  3.40  m^/sec  and  2.19  m^/sec^,  respectively,  for  the 
original  spectrum. 

The  significance  of  representing  the  shape  of  a given  wave  spectrum 
by  the  sum  of  two  symmetric  spectra  is  that  the  covariance  matrix  given 
in  Eq . (3.5)  can  be  drastically  simplified  for  each  spectrum  such  that 
all  covariances  between  xc  and  xg  become  zero  (see  Appendix  C).  That 
is,  the  covariance  matrix  of  waves  for  each  symmetric  spectrum  becomes 


/ U01 

0 

v . 
01 

“Vli 
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0 
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U2i 

vli 

V2i 
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Uoi 

0 

0 

0 

0 

0 

'Vli 

V2i 
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(3.8) 

Z-  = 
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Ol 
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0 

0 

_Vli 

V2i 

0 

y2i  i 

where  i = 1 

and 

2,  and  elements  in 

the 

above 

covariance  matrix  are 

evaluated  for  each  decomposed 

spectrum. 

As 

can 

be  seen  in  the 

covariance  matrix, 

t 

the 

• 

sets 

of 

random 

variables 

(xcl’  Xcl’  Xc2’  xc2) 

and  (x  , x , 
si  si 

’ Xs2' 

Xs2 

) are 

uncorrelated 

and 

this  makes  further 

analytical  development  of  the  distribution  function  feasible. 

We  now  consider  two  sets  of  random  variables  (x  , , x . x x ) 

cl’  cl’  c2’  c2 

i » 

and  (xsi>  xsi>  xs2*  Xs2^  ^oth  normally  distributed  with  zero  mean  and 
assumed  to  be  statistically  independent.  Hence,  the  sum  of  these  two 
sets  of  random  variables  has  a normal  distribution  with  zero  mean  and  a 
covariance  matrix  which  is  equal  to  the  sum  of  the  two  covariance 
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matrices  evaluated  for  the  symmetric  spectra.  Hereafter,  each  element 

of  the  covariance  matrix  given  in  Eq.  (3.5)  represents  the  sum  of  the 

two  elements  in  Eq . (3.8)  computed  for  each  symmetric  spectrum.  All  q in 

the  covariance  matrix  given  in  Eq.  (3.5)  are  zero. 

Since  the  inverse  of  the  covariance  matrix  can  be  evaluated,  the 

joint  probability  density  function  of  (x  , x , x ) can  now 

C -L  c i s 2 

be  written  as  follows: 


K 

2 


(3.9) 


(2tt)4  W 


where 


00  < all  x,  x < 00 
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+ 2M  ' 


, ,(.x  ,x  . 
13  cl  c2 
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(3.12) 
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The  detailed  derivation  of  Eq.  (3.9)  is  given  in  Appendix  D. 
Next,  by  applying  the  polar  coordinates 


x = R 

cl  1 

cos 

el 

xsl  ■ E1 

sin 

91 

xc2  ' E2 

cos 

92 

(3.13) 

xs2  * E2 

sin 

02‘ 

the  random  variables  in  Eq . 

(3.9)  are  transformed  to  a set 

of 

new  random 

» 

variables  (R^  R R 

V 

®1>  Qp  02,  ®2).  Since  R^  and 

r2 

represent 

the  envelope  of  a random  process  x(t)  at  times  t1  and  t2,  respectively, 
we  set  Rx  = R2  = a for  the  problem  at  issue.  Then,  the  joint  proba- 
bility density  function  of  R^,  R^,  0^ , 0^,  0^,  and  0 2 for  a specified  a 
can  be  written  by 


£<V  V V V S2>  V a)  ■ 


1*1 


L 

2lAf 


(3.14) 


0 i V "»  " ” < R2  < 0,  0 < 0j,  92  < 2i,  - » < 0 0 < 


where 


2 

L = 2oc  (Mn+  Mu  cos  (0^  02)}  + 2cc  (M^-  M^cos  (0^  02)}  (Rf  R£) 


+ M22(rJ+  R2)  + 2M24R1R2  cos(01-  02)  + M22a2(0^+  02) 
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+ 2M24ot2  0^2  cosCGj-  92)  + 2M24a  (R^-  R2  0^  sinO^  02) 

+ 2Mua2  (01  + 02)  sinCOj-  02)  (3.15) 


It  is  noted  that  Mj^,  ^2’  etc.  in  Eq . (3.15)  are  functions  of 
Vo*  V1  an<^  v2  as  shown  in  Eq.  (3.12),  which  in  turn  are  functions  of  the 
time  duration  associated  with  the  envelope  exceeding  the  level  a. 

Next,  Eq.  (3.14)  is  integrated  with  respect  to  and  0 . The 
exponential  part  of  Eq.  (3.14)  consists  of  many  terms  containing 
cos ( 0J-  e2)  and  sin  ^ - 02)  as  shown  in  Eq.  (3.15);  hence,  the 
integration  is  carried  out  with  respect  to  0^  and  (0^-  0 ).  The  results 
yield 


where 


f (R ^ > ^2  * > 


02;  a)  - 


4 

a 


( 2 TT  ) 3 | A| 


(3.16) 


P = Ifip  j_2  (M11+  M13  cos  ♦>  a2+  2^Mi2"  Ml4  cos  *>  a (Rj-  V 
+ M22(Rj  + R2)  + 2M24  R.^2  cos  <p  + M22a2(02+  02) 

+ 2M24a2  cos  <p  + 2M24«  sin  * (R^-  R^) 

+ 2M14a2  (0L+  02)  sin  <j>l  (3.17) 


After  much  lengthy  mathematical  manipulation,  the  details  of  which 
are  given  in  Appendix  E,  the  joint  probability  density  function  f(R^, 
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» » t f 

R2’  0r  02;  is  integrated  with  respect  to  Q^nd  0 , and  then,  the 
joint  probability  density  function  f^,  a)  can  be  written  as 


f(V  V “>  ' <2f>2  / — ~~  • [-  [ — ^ 


COS  (j>) 


2 2 
y - v 


+ ~ °2n[  (R1+  2 p ^i^2+  “ 2kO-P)(R,-  R»-  k))(d  4>  (3.18) 

2M22(1-P  ) 1 z 
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cos  <j> 


(3.19) 


In  order  to  carry  out  further  integration  necessary  for  the 


derivation  of 

the 

probability  density  function 

of 

the  time 

duration,  f(r  ), 
’ a+'’ 

the 

following  transformation  of  the 

random 

variables 

is  made: 

(3.20) 


Then,  the  joint  probability  density  function  of  U and  V for  a 
specified  a becomes 


2tt 
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x exp 


a (y  -v  cos  <j>) 

H 

2 2 
u - v 
o o 


1 


2 ( 1 — p ) 


(u2-  2p  u v + v2) 


0 < u < 00 
0 < v < 00 


By  applying  the  property  associated  with  the  bi-variate  normal 
probability  distribution  (Cramer,  1946),  the  2nd  term  of  the  exponential 
expression  in  Eq.  (3.21)  can  approximately  be  written  as  follows  (see 
Appendix  F): 
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(u2- 


2(1 — P ) 


2 p u v + 


v2) 


U2+ 


(3.22) 


This  approximation  method  was  used  by  Tikhonov  (1956)  in  the 
derivation  of  the  probability  distribution  of  time  duration  associated 
with  excursion  of  Gaussian  noise.  The  joint  probability  density 
function  f(u,v;  ct)  then  approximately  becomes 
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By  transformation  of  the  random  variables  R^and  R2  to  U and  V, 
respectively,  the  numerator  of  Eq . (3.1)  can  be  expressed  in  terms  of  U 
and  V as  follows: 


N (x  ) = 
a+v  a+' 


ti  i i 


RjR2  f(RL,  R2;  a)  dRxdR2 
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(3.24) 


Then,  from  Eqs . (3.23)  and  (3.24),  we  have 
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Since  the  terms  associated  with  u and  v in  Eq . 
symmetric,  we  can  simply  write 


(3.25)  are 
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du  dv 


(3.26) 
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u exp 


By  letting 
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(3.27) 


the  integration  given  in  Eq.  (3.26)  can  be  obtained  as 
r -d-  2 Y u)  „ 2 


u e 


du  = 1 + /2tt  y e2^  ^1  + $ (/2~  y) 


where 
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Z -t2 
e dt 


(3.28) 


Hence,  Eq.  (3.25)  can  be  evaluated  by 
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(3.29) 


The  first  term  of  the  above  equation  represents  the  expected  number 
(per  unit  time)  of  the  envelope  downward  crossing  with  velocity  R^ 
followed  by  the  upward  crossing  with  velocity  R„,  namely  N (t  ) 
(see  Appendix  G) . Since  N (t  ) is  smaller  than  N (t  ,)  for  a 
large  a,  the  first  term  of  Eq.  (3.29)  can  be  neglected.  Thus,  we  can 
write 


On  the  other  hand,  the  denominator  of  Eq.  (3.1)  becomes 


(3.31) 


o 

As  stated  in  regard  to  Eq.  (3.1),  the  equation  is  equivalent  to  the 
probability  density  function  of  the  time  duration,  t^+,  associated  with 
the  envelope  exceeding  a specified  level  a.  Thus,  from  Eqs.  (3.30) 
and  (3.31),  the  probability  density  function  of  ia+  can  be  obtained  as 
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(3.32) 


It  is  assumed  that  $ is  relatively  small  for  large  a , and  thereby 

2 

we  have  the  approximation  cos  n 1 - $ /2.  Then,  Eq.  (3.32)  may 

be  written  as 
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» 

x Y 


/~T  Y $(/T  Y ) 


(3.33) 


where 
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normalization  factor  to  make  the  area  under  the 
function  unity 

OO 
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(3.34) 

density 
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3.2  Probability  Distribution  of  Time  Duration 
Associated  With  Wave  Group 

In  this  section  the  probability  density  function  of  the  time 

duration  associated  with  wave  groups  will  be  derived.  The  probability 

density  function  derived  in  the  previous  section  is  that  applicable  to 

the  time  duration  that  the  envelope  of  a narrow-band  random  process 

exceeds  a specified  level.  If  the  time  duration  is  not  sufficiently 

then  there  may  be  no  wave  crest  or  only  one  wave  crest  occurring 

during  the  time  the  envelope  exceeds  the  specified  level,  and  this 

cannot  be  considered  as  a wave  group.  Therefore,  the  probability 

density  function  of  the  time  interval  associated  with  the  envelope 

exceeding  a specified  level  should  be  modified  so  that  it  is  only 

concerned  with  the  time  interval  during  which  two  or  more  wave  crests 

are  present.  For  this,  let  us  first  examine  the  relationship  between 

the  time  duration  x and  the  number  of  wave  crests  involved  in  x . 

a+ 

From  the  assumption  of  a narrow— band  random  process,  the  time 

interval  between  two  successive  wave  crests  associated  with  wave  groups 

is  considered  to  be  equal  to  the  average  zero-crossing  period,  T which 

o 

can  be  evaluated  by  the  following  formula: 


= 2 xr  / m T 
o 


m0 


(3.35) 


where 


raj  - j— th  moment  of  wave  spectrum. 
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Then,  we  can  evaluate  the  number  of  waves  during  the  time 
interval  xa+  based  on  the  average  zero-crossing  period  as  follows: 

tct+  i T0:  the  time  duration  of  the  envelope,  xa+,  is  less 

than  T , there  is  either  no  wave  crest  or  only  one  wave  crest  in  x 

a+ 

Hence,  this  situation  cannot  be  considered  as  a wave  group  although  the 
envelope  exceeds  the  specified  level  a.  The  probability  of  occurrence 
of  this  situation  can  be  evaluated  from  Eq.  (3.34)  as 


Po  = 


l 

/' 


f(T  J 
a+' 


dx 


a+ 


(3.36) 


(li)  To  < xa+  < 2 T : If  the  time  duration  x is  between  T 

o 

and  2 Tq,  then  there  is  either  one  wave  crest  or  two  wave  crests  during 

Ta+*  Tlle  probabilities  of  one  and  two  wave  crests  during  xa+  can  be 
evaluated  by 


Pr 


One  wave  crest 


Pr  ] Two  wave  crests 


(3.37) 


(3.38) 


Since  one  wave  crest  in  xa+  cannot  be  considered  as  a wave  group, 

this  situation  should  be  eliminated  from  the  probability  density 

function  for  the  time  duration  associated  with  wave  groups.  The 

probability  of  occurrence  of  one  wave  crest  in  T < x < 2 T can  hn 

o = a+  = o 

evaluated  by 


■■-/’’I1- 1 H)j 


f<v> 


(3.39) 
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(Hi)  ^ Tq  < Ta+  £ 3 Tq:  Under  this  condition,  there  are  either 
two  or  three  wave  crests  in  time  duration  Ta+;  hence,  this  situation 
can  be  considered  a wave  group. 

Thus,  in  summary,  an  exceedance  of  the  envelope  above  the  specified 
level  does  not  necessarily  mean  wave  groups,  since  there  may  be  no  wave 
crest  or  only  one  crest  if  the  exceeding  time  duration  is  not 
sufficiently  long.  The  probability  of  occurrence  of  this  situation  is 
given  by  (pQ  + p^)  as  shown  in  the  above.  Therefore,  the  probability  of 
wave  group  when  envelope  exceeds  the  specified  level  a is  given  by  the 
following  formula: 


Pr  | Wave  groups  when  envelope 
/ exceeds  a level  a 


J 


- = 1 - 


(po+  pl> 


= 1 


T 

■ir 


2T 


"V  dtc+  + 
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-1  ' 2 l1  " vjf  f(Ta+)  dTa+ 


(3.40) 


The  probability  density  function  of  time  duration  associated  with 

wave  groups  can  be  obtained  from  the  probability  density 

function  f(Ta+)  by  eliminating  the  portion  in  which  either  no  wave  crest 

or  one  crest  occurs  in  Ta+.  That  is,  the  entire  probability  density 

function  for  Ta+  £ Tq  is  discarded,  and  the  magnitude  of  the  density 

function  for  Tq  < Ta+  < 2Tq  is  modified  by  multiplying  by  the 

probability  given  in  Eq . (3.38).  This  is  shown  in  the  pictorial  sketch 

given  in  Figure  3.3.  Then,  the  modified  probability  density  function  is 

normalized  so  that  the  area  under  the  density  function  becomes  unity. 

The  probability  density  function  thus  derived,  denoted  by  f(x  ),  which 

G 
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Figure  3.3  Sketch  showing  the  relationship  between  probability  density 

function  associated  with  wave  envelope,  f (t  ),  and  that 

associated  with  wave  groups,  f(  t ).  a+ 

G 
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is  applicable  to  the  time  duration  associated  with  wave  groups  is 
summarized  as  follows: 


/ 


f(0  =< 
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l - 


(p0+  Pl)  2 l1 


f(V> 


Ta+  TG 


for  0 < t < T 
= G = o 
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for  2 T < x (3.41) 

O — Vj 


The  expected  (average)  time  duration  associated  with  wave  groups, 
denoted  by  tG  , can  thereby  be  evaluated  as 
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(3.42) 
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3.3  Number  of  Wave  Crests  In  a Group 


Once  the  probability  function  of  the  time  duration  associated  with 
wave  groups  is  derived,  statistical  information  on  the  number  of  waves 
involved  in  a group  can  be  obtained  therefrom.  For  example,  the 
probabilty  of  occurrence  of  a specified  number  of  waves  in  a group  can 
be  evaluated  from  the  probability  density  function  developed  in  the 
previous  section. 

Let  us  first  evaluate  the  average  (mean)  number  of  wave  crests  in  a 
wave  group  for  a given  wave  spectrum.  Assuming  that  the  time  interval 
between  two  successive  wave  crests  is  equal  to  the  average  zero-crossing 
period,  the  average  number  of  wave  crests  can  be  obtained  from 
Eqs.  (3.35)  and  (3.42),  as 


Next,  let  us  evaluate  the  probability  of  two  wave  crests,  three 
wave  crests,  etc.,  in  a group.  In  general,  for  the  time  duration 


(m-1)  waves  or  m-waves  in  tg  . The  probability  of  m-waves  in  a 
specified  is  given  by 


Average  number  of  wave  crests  = t / T 


(3.43) 


(m-1)  Tq  < tg  < mTQ,  where  m is  an  integer  2,3, , there  are  either 


where 


(m-1)  T < t < mT 

O — \j  — 


o 


Tq  = average  zerc 


average  zero-crossing  wave  period 
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The  probability  of  (m-1)  waves  for  a given  x becomes 

G 


Pr  < (m-1 )-waves 


M ' 1 l1 


(3.45) 


From  the  conditional  probability  given  above,  the  probability  of 

occurrence  of  m-waves  in  x , where  (m-l)T  < x„  < mT  , can  be 

« o = G = o 

evaluated  by 
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t T \ 

< m-waves 
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'-t;) 

i 
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(3.46) 


(m-l)T 


Similarly,  we  have 


mT 


Pr < (m-1)  waves>= 


1 - " l1  -~c)  ff<V  dTG  (3'47) 


(m-l)T 


By  applying  the  formulae  derived  in  Eqs.  (3.46)  and  (3.47),  we  can 
evaluate  from  the  wave  spectrum  the  probability  of  occurrence  of  a 
specified  number  of  waves  in  a group  as  follows: 


f(TG)dTG 
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(3.48) 


3.4  Example  of  Application 


As  an  example  of  the  application  of  the  prediction  method  for 
evaluating  the  time  duration  and  number  of  waves  associated  with  wave 
groups,  numerical  computations  were  carried  out  using  the  two-parameter 
wave  spectrum  shown  in  Figure  3.2  having  a significant  wave  height  of 
10.0  m and  a modal  frequency  of  0.42  rps. 

The  probability  density  function  of  the  time  duration  of  the 
envelope  exceeding  levels  of  5.0,  6.0,  7.0,  and  8.0  m above  the  mean 
water  level  was  computed  by  Eq.  (3.33),  and  the  results  are  shown  in 
Figure  3.4.  The  average  values  of  the  time  duration  x were  then 


and  the  results  compared  with  those  computed  to  date  by  the  following 


evaluated  from  Eq.  (3.33)  by 


(3.49) 


o 


formula: 


(3.50) 


PROBABILITY  DENSITY  IN  l/SEC 
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Figure  3.4  Probability  density  function  of  time  duration  associated 
with  wave  envelope  at  various  levels. 
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where  f(a,&)  is  the  joint  probability  density  function  of  amplitude  and 
velocity  of  the  envelope,  and  Uq  and  u2  are  as  defined  in  Eq.  (3.6). 

As  can  be  seen  in  Table  3.1,  the  average  time  duration  for  a given 
level  evaluated  by  the  probability  density  function  derived  in  the 
present  study  (Eq.  3.49)  is  substantially  greater  than  that  computed  by 
the  formula  given  in  Eq . (3.50). 

Additional  computations  were  carried  out  using  wave  data  measured 
at  sea.  The  data  were  taken  in  North  Sea  off  Norway  by  the  Norwegian 
Marine  Technology  Research  Institute.  The  water  depth  at  the  measured 
site  is  230  meters.  Two  sets  of  data  (significant  wave  heights  of  8.13 
m and  7.44  m)  were  analyzed.  Figure  3.5  shows  the  measured  spectrum  as 
well  as  the  sum  of  the  two  symmet  ric— shaped  spectra  used  in  the 
computations  for  each  record. 

Since  the  observation  time  for  each  record  for  which  the  seas  were 
considered  to  be  steady- state  is  not  sufficiently  long,  the  sample  size 
of  the  envelope  exceeding  a specified  level  is  rather  small.  Therefore, 
the  average  values  obtained  from  the  measured  data  may  not  be 
accurate.  Nevertheless,  a consistent  trend  can  be  seen  in  Table  3.2 
where  the  measured  and  computed  average  time  durations  above  various 
levels  are  tabulated.  That  is,  although  the  average  values  obtained 
from  the  records  are  greater  than  those  computed  by  Eq.  (3.49),  the  two 
are  very  close  for  high  crossing  levels.  On  the  other  hand,  the 

computed  values  by  Eq.  (3.50)  are  smaller  than  the  values  obtained  from 
the  records  by  a considerable  amount. 

As  stated  in  connection  with  Eq.  (3.40),  an  exceedance  of  the 
envelope  above  a specified  level  does  not  necessarily  mean  a wave 
group.  The  probability  of  occurrence  of  wave  when  the  envelope  exceeds 
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Table  3.1 


Comparison  of  average  time  duration  associated 
with  the  envelope  exceeding  a specified  level. 


Level 

Average  time  duration  computed  by 

Eq.  (3.49) 

Eq.  (3.50) 

5.0  m. 

8.2  sec. 

5.6  sec. 

6.0 

7.3 

4.6 

7.0 

6.4 

4.0 

8.0 

5.5 

3.5 

SPECTRAL  DENSITY  IN  M2-SEC  SPECTRAL  DENSITY  IN  M2  -SEC 
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(a)  Wave  spectrum  A,  significant  wave  height  8.13  m. 


(b)  Wave  spectrum  B,  significant  wave  height  7.44  m. 
Figure  3.5  Measured  wave  spectra  used  in  analysis. 


Table  3.2  Comparison  between  observed  and  computed  average  time  durations 
associated  with  the  envelope  exceeding  a specified  level. 


Level 

Observed 

Computed 

average  time  duration 

Sample 

size 

Average 
time  duration 

Eq. (3.49) 

Eq. (3.50) 

< e' 
co 

E t-h 

4.0  m. 

9 

15.9  sec. 

9.1  sec. 

3.9  sec. 

3 • 

u oo 

4J 

CJ 

0)  II 
a 

CO  w 

Sd 

5.0 

6 

10.8 

7.9 

3.1 

6.0 

2 

7.0 

6.7 

2.6 

m e 

4.0  m. 

7 

10.7  sec. 

7.4  sec. 

3.7  sec. 

trum 

7.44 

5.0 

4 

7.1 

6.3 

3.0 

a 

a)  ii 
a w 
CO  jc 
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group.  The  probability  of  occurrence  of  wave  when  the  envelope  exceeds 
the  specified  level  was  calculated  by  Eq.  (3.40)  using  the  two-parameter 
wave  spectrum  and  the  results  are  tabulated  in  Table  3.3. 

As  can  be  seen  in  the  table,  the  majority  of  the  level  crossings 
cannot  be  considered  as  indicative  of  the  occurrence  of  wave  groups 
since  there  is  either  no  wave  crest  or  only  one  crest  during  the 
crossing  time  interval. 

Figure  3.6  shows  the  probability  density  function  of  the  time 
duration  associated  with  wave  groups  exceeding  levels  of  5.0,  6.0,  7.0, 
and  8.0  m above  the  mean  water  level  computed  by  Eq.  (3.41)  using  the 
two-parameter  wave  spectrum.  It  can  be  seen  from  a comparison  of 

Figures  3.4  and  3.6  that  there  is  a substantial  difference  in  these 
probability  density  functions. 

Table  3.4  shows  the  comparison  between  observed  and  computed 
average  time  durations  associated  with  wave  group  for  the  measurerd  wave 

spectrum  A shown  in  Figure  3.5(a).  Agreement  between  them  appears  to  be 
satisfactory. 
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Table  3.3 


Probability  of  occurrence  of  wave  group  when 
envelope  exceeds  a specified  level. 


Probability 

of 

No  wave  group 

Level 

Wave  group 

No  wave  crest 

One  wave  crest 

5.0  m. 

0.16 

0.71 

0.13 

6.0 

0.12 

0.76 

0.12 

7.0 

0.08 

0.83 

0.09 

8.0 

0.05 

0.88 

0.07 

PROBABILITY  DENSITY  IN  I /SEC 
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Figure  3.6  Probability  density  function  of  time  duration  associated 
with  wave  groups  at  various  levels. 


Table  3.4  Comparison  between  observed  and  computed  average 
time  durations  associated  with  wave  group 
(Spectrum  A,  H = 8.13  m.) 


Level 

Observed 

Computed  average 
time  duration 

Sample 

size 

Average 
time  duration 

4.0  m. 

6 

19.0  sec. 

17.9  sec. 

5.0 

1 

13.7 

17.0 

CHAPTER  IV 

FREQUENCY  OF  OCCURRENCE  OF  WAVE  GROUPS 


4,1  -Probability  Distribution  of  Time  Interval  Between  Successive 
Crossings  of  Specified  Level  by  Envelope 

We  first  consider  the  level  crossing  of  the  envelope  of  a narrow 
band  random  process  and  derive  the  probability  distribution  of  the  time 
interval  between  successive  positive  crossings  of  the  envelope,  R(t),  at 
a specified  level  a which  is  denoted  by  xa  as  shown  in  Figure  1.3. 

The  distribution  function  of  xa  can  be  evaluated  by  applying  the 
concept  presented  in  Chapter  III  of  the  present  study.  That  is,  the 
probability  distribution  function  of  time  duration  of  the  envelope 
exceeding  a specified  level  a,  denoted  by  xa+,  was  developed  in 
Chapter  III.  In  a similar  fashion,  the  probability  distribution 
function  of  time  duration  of  the  envelope  below  a specified  level, 
denoted  by  xa_,  may  be  developed.  Since  the  time  interval  between 
successive  positive  crossings  of  the  envelope,  x is  the  sum  of  x . and 

CX  T* 

ta_,  the  probability  distribution  function  of  xa  can  be  derived  by  the 
convolution  integral  of  two  probability  density  functions  f(xa+)  and 
f(Ta_).  This  approach,  however,  has  a drawback  in  the  accuracy  of  the 
probability  density  function  f(xa_).  To  elaborate,  the  following 
explanation  is  given: 

The  method  for  deriving  the  probability  density  function  developed 
in  Chapter  III  is  valid  for  a relatively  small  time  duration.  Since  a 
fairly  high  level  of  envelope  crossings  is  considered  for  wave  groups, 
the  time  duration  xa+  is  small,  but  the  time  duration  of  the  envelope 
crossing  below  the  level,  xa_,  is  essentially  large.  Hence,  there  is  no 
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assurance  as  to  the  accuracy  of  the  probability  density  function  f(TQ_) 
if  it  is  derived  following  the  approach  developed  in  Chapter  III. 

One  way  to  overcome  the  difficulty  involved  in  the  derivation  of 
the  probability  density  function  of  the  time  interval  between  successive 
positive  crossings  of  the  envelope  is  to  assume  that  the  probability 
density  function  is  approximately  equal  to  that  of  the  time  interval 
between  successive  maxima  of  the  envelope  above  a specified  level  a, 
denoted  by  Tm  in  Figure  1.3.  This  assumption  is  permissible  because  the 
crossings  take  place  at  a fairly  large  distance  above  the  zero-line  for 
wave  groups,  and  this  results  in  the  time  interval  between  positive 
crossings  being  nearly  equal  to  that  between  successive  maxima  of  the 


envelope. 


Following  Rice's  method  developed  for  evaluating  the  avei 
of  maxima  per  unit  time  (Rice,  1945),  we  assume  waves  to  be 
band  random  process  whose  profile  x(t)  can  be  written  as 


average  number 


a narrow- 


x(t)  = Z c 
i n 


(4.1) 


where 


x (t)  - Z c cos{(oj  - m)t  - e } 
r.  = i n n nJ 


xs(t)  - Z sin{(wn  - w)t  - e } 


(4.2) 


c 


— J 2 S ( a)  ) \ Am 
n i 


n 


mean  frequency  of  wave  spectrum 
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Since  we  are  interested  in  the  maxima  of  the  envelope  shown  in 
Figure  1.3,  we  consider  a random  vector  X which  is  composed  of  a set  of 
six  random  variables  xc>  xc,  xc>  xg)  xg,  xg  representing  the  cosine 
and  sine  components  of  displacement,  velocity,  and  acceleration. 

Each  element  of  X obeys  the  normal  probability  distribution  and  we 
can  write  the  probability  density  function  of  X in  the  following  form: 


f(X)  = 


1 


1 


- j X’lf 1 X 


( 2 TT  ) 


1/2 


(4.3) 


where 


X'  " (V  V V v v V 


E is  the  covariance  matrix  given  by 


E = 


/"< 

/ 0 

-y, 

0 

\: 


0 

U2 

0 

0 

0 

V- 


~\lr 

0 

0 

‘U3 

0 


0 

0 

0 

y 

c 

0 

-U0 


0 

0 

0 

hn 


(4.4) 


where 


CO 

v-I 


(w  - w) r . S(uj)  da) 


(4.5) 


S(u>)  = wave  spectrum 


The  derivation  of  elements  of  the  covariance  matrix  associated  with 
displacement  and  velocity  was  given  in  Appendix  A.  The  derivation  of 
the  elements  associated  with  acceleration  is  given  in  Appendix  H. 

As  can  be  seen  in  Eq.  (4.4),  it  is  difficult  to  invert  the  covari- 
ance matrix  E as  needed  in  Eq . (4.3). 


However,  by  interchanging  the 
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positions  of  x and  x the  matrix  given  in  Eq.  (4.4)  can  be  written 


as 


E = 


0 

-M, 

0 

0 

0 


0 

P, 

-p, 

0 

0 

0 


~W2 

U4 

0 

0 

0 


0 

0 

0 

p 

o 

0 

-p„ 


0 

0 

0 

0 

P„ 


0 
0 
0 

-p 


(4.6) 


and  thereby  the  inversion  of  the  matrix  is  feasible. 

The  joint  probability  density  function  of  (xc,  xg , x£,  xg,  x , x ) 


can  now  be  written  as  follows: 


T " f II 

f(XC,  Xg,  Xc,  Xg,  Xc,  Xg)  = 


l 


K 

1 2 

TiTF  JaT  6 

I II 

-co  < all  x,  x,  x < 00 


(4.7) 


where 


|A|  - „o(w2  - II2  ) - 


K 


_ 1 ru  , 2 2.  , ' 2 '2  "7  " 

' T^T  ^ 11(XC  + X«)  + + O + Moq(x2  + X 


22  c 


33  c 


+ 2M12<xc  xs  - xc  xs)  + 2M13(xc  xc  + xs  xs>  + 2M23(xc  xs  ‘ xc  xs>) 


(4.8) 
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and 


M11  = W2  " ^3 


M22  = % M4  ' U2 


M33  Uo  U2 


M12  U2  U3 


M13  U2 


= U M- 
23  o 3 


(4.9) 


The  detailed  derivation  of  Eq . (4.7)  is  given  in  Appendix  I. 

Next,  by  applying  the  following  polar  coordinates,  the  joint  proba- 
bility density  function  given  in  Eq.  (4.7)  will  be  transformed  to  that 

in  in 

of  R,  R,  R,  0,  0,  0.  That  is,  by  letting 


xc 

— 

R 

COS 

0 

xs 

= 

R 

sin 

0 

1 

xc 

= 

I 

R 

cos 

0 

- 

R 

0 

sin  0 

1 

xs 

= 

t 

R 

sin 

0 

+ 

R 

§ 

cos  0 

If 

xc 

= 

ii 

R 

cos 

0 

- 

2 

I 

R 

1 

9 sin 

If 

xs 

= 

ii 

R 

sin 

0 

+ 

2 

I 

R 

t 

6 cos 

9 - R 9^  cos  0 - R 0 sin  0 
0 - R 0^  sin  0 + R 0 cos  0, 


(4.10) 


the  joint  probability  density  function  can  be  written  by 


_ , I II  | || 

f(R,  R,  R,  9,  0 , 0)  = 


L 

RI * 3  1 

(T^FIaT  e 


0 < R < ”,  < R < «, 


I I 

-°°  < R < 0,  0 < 9 < 2tt 

I If 

-C°  < 0 < 00,  _oo  < 0 < oo> 


(4.11) 
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.2  . ..  2 . „ 2 12,  . ..  ”2  ,'2  12  2 '4 


and  L = MnR  + M^R  + R 0*)  + M3;J(RZ  + 4RZ  0Z  + R2  0 


2 "2  ii  to  t t ti  o , „ 

+ R 0 - 2 RR  0 + 4RR  0 0)  + 2 M12  RZ  0 + 2 M13(RR 


2 ' 2 
R 0 ) 


tt  t ' o i o t o t tt 

+2M23(RR  0 - 2RZ  0 - RZ  0J  - RR  0) 


(4.12) 


By  integrating  Eq.  (4.11)  with  respect  to  0 and  0,  we  have 


f(R,  R,  R,  0)  = 


(2ir)3/2(  | a|  M33)1/2 


P 

2|aT 


(4.13) 


0 * R < 00 , -°°  < R < °°,  -«  < R < 0,  -oo  < e < 


where 


P - MnR2  + M22(R2  + R2  02 ) + M33(R2  + 4 R2  02  + R2  04  - 2 R R 02) 


+ 2M12R  0 + 2 Mn(R  R - R2  02)  + 2M23(R  R 0 - 2 R2  0 - R2  03) 


+ (R2/M33)(2  M3;j  0 + M23)2. 


(4.14) 


Next,  the  probability  density  function  of  positive  maxima  of  the 
envelope  and  phase  velocity  will  be  derived  from  Eq . (4.13).  Let  E be  a 
random  variable  representing  the  local  positive  maxima  of  the  envelope. 
Then,  the  probability  that  E will  exceed  a level  ? with  a phase  velocity 
& is  equivalent  to  the  expected  value  of  the  ratio  of  the  number  of 
positive  maxima  above  the  level  5 per  unit  time  with  velocity  9,  denoted 

by  N(g,$),  to  the  total  number  of  positive  maxima  per  unit  time,  denoted 
by  N(5+).  That  is, 
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t 

Pr{H  > ?,  0}  = 1 - F(5,  0)  = E[^->9)]  (4. 15) 

where  E[  ] stands  for  the  expected  value  and  F(5,0)  is  the  joint  cumula- 
tive distribution  function  of  £ and  0. 

If  we  assume  that  N(5,0)/N(g+)  and  N(?+)  are  statistically  indepen- 
dent, then  Eq.  (4.15)  can  be  written  as 


Pr{=  > 5,  8}  = 1 - FCC,  9) 


(4.16) 


Taking  R 0 at  positive  maxima  into  consideration,  the  numerator 
and  denominator  of  Eq.  (4.16)  can  be  evaluated  as  follows: 


e oo  o 

E[N(5,  0)]  = N^'  = -f  f J R f(R.  0,  R,  0)  dR  dR  d0  (4.17) 


E[N(?+)]  = N 


From  Eqs.  (4.16),  (4.17),  and  (4.18),  the  joint  probability  density 
function  of  the  envelope  maxima  and  phase  velocity  0 can  be  obtained  as 


00  00  V 

■■/// 

— OO  0 _ao 


R f (R,  0,  R,  0)  dR  dR  d0 


(4.18) 


where  0 < E,  < 


«—00 


< b < «. 
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The  numerator  of  Eq.  (4.19)  can  be  evaluated  analytically;  however 
the  integration  with  respect  to  R and  0 involved  in  the  denominator  can 
be  evaluated  numerically.  The  details  of  the  derivation  are  given  in 

Appendix  J.  The  resulting  joint  probability  density  function  of  § and  0 
becomes 


f(5,  9)  = 


2 2 
?2  e ^ u[l  + 5 v/tt  {l  + $(£V)}] 


r oo  2 ~ , 

J^J  Q R e R U[1  + R v/tt  e^Rv^  {l  + $(RV)}]  dR  d0 


0 < 5 < °°,  -»  < 0 < 


(4.20) 


where 


U ' 2f Af  tMll  + 2 M12  9 + (M22  " 2 M13}  '°2  ~ 2 M23  ^ + M33 


v = 


t2l*l  M33} 


m <M13  + M23  6 - M33  6 > 


(4.21) 


From  the  joint  probability  density  function  of  5 and  0,  we  can 

derive  the  conditional  probability  density  function  of  0 given  £.  That 
is, 


f(0|?) 


f(5,  9) 


f(£,  0)  d0 


(4.22) 
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Next,  the  conditional  probability  density  function  given  in 
Eq.  (4.22)  will  be  converted  to  the  conditional  probability  density 
function  of  the  time  interval  between  successive  peaks  above  a specified 
level  a by  letting  E,  = a and 

6 = 2Trk/ra,  (4.23) 

where  the  parameter  k is  determined  by  the  iteration  method  from  a con- 
sideration of  the  probability  of  envelope  crossing  at  the  level  a.  The 
derivation  of  Eq.  (4.23)  is  given  in  Appendix  K along  with  a method  to 
determine  the  parameter  k.  Thus,  from  Eqs.  (4.20)  through  (4.23),  the 
conditional  probability  density  function  of  xa,  given  E,  = a,  can  be 
written  as  follows: 


f(TJa) 


e~^u  [1  + 5 v/i  eav)2 
e"?2u[l  + § v/i  e(?v)2 


{1  + *(5v)}] 

{l  + <KCv)}]  d9 


2 irk 


-a2u' 

e 


[l  + av'/if  e^av  ^ {l  + $ ( av ' ) } ] 
+ cxv/tt  e^av^  {l  + $(av)}]  d9 


C = a 

? 

9 = 2trk/T 

a 


2irk 


2 ’ 
T 

a 


(4.24) 


where 


(4.25) 
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constant  for  normalization  of  the  probability  density  function 


We  may  recall  the  assumption  that  the  time  interval  between  suc- 
cessive peaks  of  the  envelope  above  a specified  level  a is  equal  to  the 
time  interval  between  successive  up-crossings  of  the  level  by  the  enve- 
lope. Hence,  the  conditional  probability  density  function  given  in  Eq. 
(4.24)  represents,  as  well,  the  conditional  probability  density  function 
of  the  time  interval  between  successive  up-crossings  of  a specified 
level  a by  the  envelope. 

Equation  (4.24)  can  also  be  expressed  in  terms  of  the  central 
moments  defined  in  Eq.  (4.5)  with  the  aid  of  Eq.  (4.9).  In  this  case, 
the  parameters  u,  u' , v and  v’  become 


u 


'5  u, 

- 2u  U-.  9 + p u 0 } 

o 3 o 2 ‘ 


* <Uo  "4  - 3 A) 


T 

a 


2 


(4.26) 
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Numerical  computations  are  carried  out  for  the  two-parameter  wave 
spectrum  with  significant  wave  height  of  10.0  m and  a modal  frequency  of 
0.42  rps  shown  in  Figure  4.1.  The  probability  density  functions  of  the 
time  interval  between  two  successive  envelope  up-crossings  for  level  of 
5.0,  6.0,  and  7.0  m above  the  still  water  level  are  computed  by 
Eq.  (4.24)  and  the  results  are  shown  in  Figure  4.2. 

4 • 2 Probability  Distribution  of  Time  Interval  Between  Wave  Groups 
As  discussed  in  the  Introduction  referring  to  Figure  1.3,  an  up- 
crossing  of  the  envelope  does  not  necessarily  mean  the  occurrence  of  a 
wave  group.  The  probability  of  occurrence  of  wave  groups  when  the 
envelope  exceeds  a specified  level  a was  obtained  in  Eq.  (3.40)  of 
Chapter  III.  Let  this  probability  be  denoted  by  p.  Then,  it  can  be 
said  that  p-fraction  of  f(ra|a)  derived  in  Eq.  (4.24)  represents  the 
probability  density  function  of  the  time  interval  between  two  wave 
groups,  but  ( l-p)-f raction  of  f(xa|a)  should  be  modified  so  that  it  is 
associated  with  the  time  interval  between  two  wave  groups.  The 
modification  of  the  probability  density  function  may  be  made  as  follows: 
First,  let  us  consider  three  up-crossings  of  the  envelope.  If  the 
first  and  third  up-crossings  are  associated  with  wave  groups  but  the 
time  duration  of  the  second  up-crossing  is  so  short  that  it  does  not 
constitute  a wave  group,  then  the  time  duration  between  two  wave  groups 
is  that  between  the  first  and  third  up-crossings.  This  is  the  situation 
of  the  up— crossings  A,  B,  and  C shown  in  Figure  1.3.  In  this  case  the 
time  interval  between  Points  A and  C is  considered  as  the  sum  of  two 
independent  random  variables,  each  having  the  probability  density  func- 
tion derived  in  Eq.  (4.24).  By  writing  the  sum  of  two  independent  time 
intervals  as  t2a,  its  probability  density  function  for  a specified  level 
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Figure  4.1  Two  parameter  wave  spectrum  (significant  wave  height 
10.0  m. , modal  frquency  0.42  r.p.s.). 
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Figure  4.2  Probability  density  function  of  time  interval  between  two 
successive  envelope  up-crossings  for  level  of  5.0,  6.0, 
and  7.0  m.  computed  for  the  two  parameter  wave  spectrum 
( significant  wave  height  10.0  m.). 
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a can  be  obtained  by  the  following  convolution  integral: 

f(Ta|a)  . f((T2a-  Ta)|cO  d Ta  (4.27) 

f(Ta|a)  = probability  density  function  of  time  interval  between 
successive  up-crossings  of  envelope  above  a specified 
level  a (see  Eq.  4.24). 


The  same  concept  can  be  applied  to  four  up— crossings  of  the  enve- 
lope in  which  the  first  and  fourth  up-crossings  are  associated  with  wave 
groups.  In  this  case  the  time  interval  between  successive  wave  groups, 
denoted  by  x3a,  can  be  considered  as  the  sum  of  three  independent  random 
variables  each  having  the  probability  density  function  given  in  Eq . 
(4.24).  Its  probability  density  function  can  be  obtained  from 
Eqs.  (4.24)  and  (4.27)  as 

£<h>>  -/3“  f(T2a 

The  same  procedure  can  be  extended  to  any  number  of  up— crossings  of 
the  envelope  in  which  the  first  and  the  last  up-crossings  are  associated 
with  wave  groups.  An  example  of  the  probability  density  function 
f(TncJa)  f°r  n from  2 to  5,  computed  for  the  two-parameter  wave  spec- 
trum, is  shown  in  Figure  4.3.  The  crossing  level  a is  5 m in  this 
example.  Included  also  in  the  figure,  for  comparison,  is  the  proba- 
bility density  function  f(xja)  shown  in  Figure  4.2. 

From  the  results  of  analysis  of  wave  records  obtained  at  sea,  it 
was  found  that,  in  practice,  an  upper  limit  of  x5a  is  sufficient  for 


a)  . f (t 


3a 


" T2a) ' 


a dr 


2a 


(4.28) 
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Figure  4.3  Probability  density  functions  f(x  |a)  for  n from  1 to  5 
computed  for  the  two  parameter  wave  spectrum  (significant 
wave  height  10.0  m. , level  a = 5.0  m.) 
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consideration  of  the  longest  time  interval  between  up-crossings  for  wave 
groups,  and  that  T3a>  T4a»  an<^  T5a  aPPear  to  have  equal  chances  of 

occurrence  for  high  level  crossings.  Here,  implies  the  time  inter- 

val between  the  first  and  the  last  up-crossings  of  a total  of  six 
crossings  of  the  envelope.  Therefore,  in  the  present  analysis,  the 
probability  density  functions  f(TnJa)  for  n from  2 to  5 are  evaluated 
through  convolution  integrals  and  then  these  probability  density  func- 
tions are  accumulated  with  equal  weight  (l-p)/4  along  with  the  proba- 
bility density  function  f(Ta|a)  which  has  a weight  p.  Thus,  in  summary, 
the  probability  density  function  of  the  time  interval  between  successive 
wave  groups  for  a specified  level  a,  denoted  by  f(raG|a),  can  be  written 
as 

f(vla)  = P • f(TJa)  + j^f(T2ala)  + f(T3cJa) 

+ f(T4cJa)  + f(T5cJa)j  (4.29) 

where  p = probability  of  occurrence  of  wave  groups  when  the  envelope 
exceeds  a specified  level  a (see  Eq.  (3.40)  in  Chapter  3). 

f<TnJ°°  ~fo  £(<V  T(n-l)a)l“)-  f ( T(„-l  )J  ° ) • dt(„-i)c, 

n = 2,  3,  4,  and  5 (4.30) 

As  an  example,  Figure  4.4  shows  the  probability  density  function 

^TaGla)  computed  for  the  two-parameter  wave  spectrum  for  the  level 
d = 5 m. 

The  average  time  interval  between  successive  wave  groups,  denoted 
^ctG»  can  be  evaluated  from  Eq.  (4.29)  as 
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Figure  4.4  Probability  density  function  of  time  interval  between 
successive  wave  groups  computed  for  the  two  parameter 
wave  spectrum  (significant  wave  height  10.0  m. , level 
oi  = 5.0  m. ) 
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x 


aG 


aG 


f(TaGl°0  dx 


aG 


(4.31) 


4,3  Evaluation  of  Occurrence  of  Wave  Groups 

It  is  of  considerable  interest  to  evaluate  for  a given  wave  spec- 
trum how  many  times  a wave  group  will  occur  in  a specified  time  period, 
say  30  minutes  or  one  hour,  as  long  as  the  sea  severity  is  unchanged. 
For  this,  renewal  theory  is  applied  to  the  probability  density  function 
of  the  time  interval  between  wave  groups. 

It  has  often  been  said  that  the  wave  group  phenomenon  obeys  the 
Poisson  random  process.  This  assumption  may  be  true  if  we  consider  the 
wave  group  phenomenon  to  be  a rare  event.  In  this  case,  the  time  inter- 
val between  successive  occurrence  of  the  event  (wave  groups  for  the 
present  problem)  must  follow  the  exponential  probability  law.  However, 
as  shown  in  Figure  4.4,  the  probability  density  function  of  the  time 
interval  between  wave  groups,  f(TaG|o),  cannot  be  considered  to  be  an 
exponential  distribution,  in  general.  Hence,  in  the  present  analysis, 
the  probability  density  function  fd^a)  is  approximated  by  the 
following  gamma  probability  distribution  which  is  a generalized  form  of 
the  exponential  distribution: 


r(m) 


,m  m-1 

x t«g 


-Ax 


aG 


0 < x < 
aG 


(4.32) 


Note  that  fdaG)  reduces  to  the  exponential  distribution  for  m = 1, 
The  parameters  ra  and  A of  the  eamma  ihscHKnH^  v-  //. 


gamma  distribution  given  in  Eq.  (4.32) 


can 
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be  determined  by  equating  the  mean  and  variance  computed  from  the  proba- 
bility density  function  given  in  Eq . (4.29)  to  the  mean  and  variance  of 
the  gamma  distribution.  An  example  of  a comparison  between  the  proba- 
bility density  functions  computed  from  Eq.  (4.29)  and  the  gamma  distri- 
bution for  the  two-parameter  wave  spectrum  is  shown  in  Figure  4.5. 

Based  on  the  probability  density  function  given  in  Eq . (4.32),  the 

probability  of  n-occurrences  of  wave  groups  in  a specified  time  period 
can  be  derived  as  follows: 

Let  Tn  be  the  waiting  time  up  to  the  n-th  wave  group.  Then,  we  may 

write 

Tn  “ (t<xG*1  + ^TotG^2  + (Ta g*3  + + <TaG>n»  (4.33) 

where  all  TaG  in  Eq.  (4.33)  are  statistically  independent  and  follow  the 
probability  density  function  given  in  Eq . (4.32).  Since  the  charac- 
teristic function  of  the  gamma  probability  distribution  given  in 
Eq . (4.32)  can  be  written  as 

v>  ■ (■  - 

the  characteristic  function  of  Tn  becomes 

Vu) '(' 

This  is  also  the  characteristic  function  of  the  gamma  probability  dis- 
tribution with  the  parameters  X and  mn.  Thus,  the  probability  density 
function  of  Tn  can  be  obtained  as 


(4.34) 


(4.35) 
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igure  4.5  Comparison  between  probability  density  function  of  time  interval 
between  successive  wave  groups  and  gamma  probability  density 
function  ( the  two-parameter  wave  spectrum  of  significant 
wave  height  10.0  m. , and  level  a = 5.0  m.). 
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f(T  ) 
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— XT 

,mn  mn-1  n 

T7  A l e 

r(mn)  n 


(4.36) 


Similarly,  the  probability  density  function  for  the  waiting  time  up 
to  the  (n+1 ) th  wave  group,  denoted  by  Tn+1,  can  be  written  as 


f(T_,)  = TT-rr—  - T Ara(n+1)  Tm(n+1  )-l  e_ATn+l 


n+1  r(m(n  + 1 ) J 


(4.37) 


Then,  we  have  the  following  relationship: 


p rn-occurrences  of  wave i , , 

groups  in  time  t J “ PrlTn  < t}  - Pr{ln+1  < t} 


= [ f (T  ) dT  - J f(T  ,) 

J0  n n 1 n+1 


dT 


n+1 


Y(mn,  At)  y(m(n  + 1),  At) 
T (mn) 


r(m(n  + 1)) 
n = 1,  2,  3,  ... 

where  y(  ) = incomplete  gamma  function. 


(4.38) 


In  particular,  for  n = 0,  we  have 


pr /no  occurrence  of  a 
wave  group  in  time 


t}  = 1 - Pr{T1  < t} 


= 1 - 


y(m,  At) 

T(m) 


(4.39) 
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As  an  example  of  application  of  the  prediction  method  derived  in 
above,  numerical  computations  are  carried  out  using  wave  data  measured 
in  the  North  Sea  off  Norway.  The  recording  time  was  17  minutes,  and  the 
significant  wave  height  was  8.13  m.  The  spectral  density  function 
obtained  from  data  is  shown  as  Wave  Spectrum  A in  Figure  3.5  in  Chap- 
ter III  of  the  present  dissertation. 

The  probability  density  function  of  the  time  interval  between  suc- 
cessive wave  group,  f(xaG)>  computed  for  a level  of  a = 4.0  m is  shown 
in  Figure  4.6  along  with  the  gamma  distribution  computed  by  Eq.  (4.32). 

The  probabilities  of  wave  groups  in  17  minutes  for  a level  of 
a = 4.0  m are  computed  for  various  number  of  occurrences  by  Eq . (4.38) 
and  the  results  are  tabulated  in  Table  4.1  As  can  be  seen  in  the  table, 
the  probability  of  occurrence  of  4,  5,  and  6 wave  groups  is  high;  speci- 
fically the  highest  possibility  is  0.204  for  5 occurrences.  This  pre- 
dicted result  agrees  very  well  with  the  measurements  in  which  five  wave 
groups  were  observed  in  17  minutes  at  the  level  of  4 meters. 

Currently  available  methods  for  evaluating  the  time  interval 
between  successive  wave  groups  consider  the  average  time  interval 
between  two  successive  envelope  crossings  of  a specified  level.  This 
results  in  the  average  number  of  occurrences  of  wave  groups  being 
extremely  large,  the  average  being  37.7  occurrences  in  17  minutes  com- 
puted for  the  example  shown  in  Table  4.1.  The  average  number  of  occur- 
rences of  wave  groups  computed  by  the  present  theory  is  5.5  in  17 
minutes.  This  number  can  be  obtained  from  the  average  time  interval 
between  successive  wave  groups  given  in  Eq.  (4.31). 
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Figure  4.6  Comparison  between  probability  density  function  of  time 

interval  between  successive  wave  groups  and  gamma  probability 
density  function  ( measured  spectrum  of  significant  wave 
height  8.13  m. , and  level  a = 4.0  m. ) 
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Table  4.1  Predicted  frequency  of  occurrence  of  wave  groups  and 
comparison  with  measured  number  of  wave  groups. 

Data  measured  in  North  Sea  off  Norway  for  17  minutes, 
significant  wave  height  8.13  m. , level  a = 4.0  m. 


Predicted  probability 
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wave 

group 
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= 0.001 
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group 
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= 0.012 
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groups 
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groups 
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groups 
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groups 

} 
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= 0.073 
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groups 
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= 0.036 

Pr  { 

more 

wave 

: than  10 
groups 
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= 0.024 

Measured  number 
of  wave  groups 


CHAPTER  V 
CONCLUSIONS 


This  study  presents  a method  for  predicting  various  stochastic 
characteristics  of  wave  groups  in  random  seas  where  the  water  depth  is 
sufficiently  deep.  The  purpose  of  this  study  is  to  develop  analyti- 
cally, for  a given  wave  spectrum,  (i)  the  probability  density  function 
of  the  time  duration  of  a wave  group,  (ii)  the  probability  of  occurrence 
of  a specified  number  of  waves  in  a wave  group,  (iii)  the  probability 
density  function  of  the  time  interval  between  successive  wave  groups  and 
(iv)  the  frequency  of  occurrence  of  wave  groups  in  a specified  time. 

It  has  been  customary,  to  date,  to  consider  the  exceedance  of  the 
wave  envelope  above  a certain  level  to  identify  a wave  group.  However, 
if  the  time  duration  above  a certain  level  is  relatively  short,  there 
may  be  only  one  wave  crest  (or  no  wave  crest)  in  the  duration,  which 
obviously  does  not  constitute  a wave  group  even  though  the  envelope 
exceeds  the  specified  level.  Hence,  in  the  present  study,  a sequence  of 
(at  least  two)  high  waves  exceeding  a specified  level  is  considered  as  a 
wave  group. 

The  probability  density  function  applicable  for  the  time  duration 
of  the  wave  envelope  above  a certain  level  is  analytically  developed 
based  on  the  concept  associated  with  the  level  crossing  problem  of  a 
narrow  band  Gaussian  random  process.  This  probability  density  function 
is  then  modified  to  that  applicable  for  the  wave  group  (Equation  3.41). 

Formula  to  evaluate  the  probability  of  the  number  of  wave  crests  in 
a group  is  also  derived  based  on  the  probability  density  function  of 
time  duration  associated  with  wave  groups  (Equation  3.48). 
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The  probability  density  function  of  the  time  interval  between 
successive  envelope  crossings  of  a specified  level  is  developed,  and 
later  this  probability  density  function  is  modified  taking  into  con- 
sideration the  probability  of  occurrence  of  wave  groups  when  the  enve- 
lope exceeds  a specified  level.  The  probability  density  function  of  the 
interval  applicable  for  the  wave  groups,  thusly  derived,  is  given  in 
Equation  4.29.  Renewal  theory  is  then  applied  to  the  probability 
density  function  of  the  time  interval  between  wave  groups  for  evaluating 
the  frequency  of  occurrence  of  wave  groups  in  a specified  time,  say  30 
minutes  or  one  hour  (Equation  4.38). 

The  average  time  duration  associated  with  the  envelope  exceeding  a 
specified  level  computed  by  the  newly  developed  theory  agrees  reasonably 
well  with  data  observed  in  the  North  Sea  for  a high  crossing  level, 
while  the  average  time  duration  computed  by  applying  the  formula  in 
current  use  is  substantially  smaller  than  observed  data  (See  Table  3.2). 
Also,  the  average  time  duration  associated  with  wave  groups  computed  by 
the  present  study  agrees  well  with  that  which  is  observed  in  the  data 
(Table  3.4). 

The  frequency  of  occurrence  of  wave  groups  predicted  by  the  present 
theory  agrees  very  well  with  that  observed,  while  the  frequency  of 
occurrence  computed  by  applying  the  formula  currently  used  is 
substantially  greater  than  observed  data  from  the  North  Sea. 


APPENDIX  A DERIVATION  OF  ELEMENTS  OF  COVARIANCE  MATRIX 

OF  THE  RANDOM  VECTOR  X'  = (x  , x , x o x x v Y * \ 

v cl*  xcl»  xc2»  XC2>  xsl*  xsl»  xs2»  xs2^ 

The  covariance  matrix  of  the  random  vector  X'  = (xc^,  xc^,  x 2>  x 2, 

» t 

xsl>  xsl>  xs2’  xs2^  can  be  written  by 
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where  x^ , xg  are  expressed  by 
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(A. A. 2) 


lx  are 


The  subscripts  i = 1,2  of  the  elements  of  the  covariance  matrix 
associated  with  the  values  of  xc(t)  and  xg(t)  at  times  tj  and  t2  • £n  are 

the  phases.  cn  are  the  amplitudes  which  are  related  with  the  spectral 
density  function  as 
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Each  element  of  the  covariance  matrix  may  be  evaluated  as  follows: 
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The  bar  denotes  the  average  with  respect  to  time  t.  By  letting  Am  small 
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Thus,  E [xclxc2]  is  equal  to  the  area  under  the  function  which  is  the 

product  of  spectral  density  function  and  cos  (to  - w)t,  and  is  denoted 
by  v . That  is, 
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Thus,  E txclxc2 ] is  equal  to  the  derivative  of  Vq  with  respect  to  t,  and  it 
is  given  in  the  integral  form  as 
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Thus,  E Uclxsl]  Is  equal  to  the  first  moment  of  the  spectrum  about  Us 
frequency  and  it  becomes 


mean 
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Thus,  E [xclxg2]  is  the  derivative  of  nQ  with  respect  to  t.  It  is  denoted 
by  n1  and  is  given  in  the  integration  form  as 
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Thus,  E [x^]  is  equal  to  the  second  moment  of  the  spectrum  about  its  mean 
frequency  which  can  be  expressed  by 
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Thus,  E f xc^ xc2 ] is  equal  to  the  derivative  of  Vq  with  respect  to  t but 
with  opposite  sign.  It  is  given  in  the  integral  form  as 
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Thus,  E [xclxc2l  is  equal  to  the  derivative  of  vx  with  respect  to  t.  It 
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Thus,  E [*clxgl]  is  equal  to  the  first  moment  of  the  spectral  density 
function  about  its  mean  frequency,  which  is  zero.  That  is, 


E [xclxsl]  = 


CO 

J 


(to  ~ to)  S(m)  dto  = y = 


(A. A. 15) 


E * - E 


Ec  (to  - 
n n 


(o)  sin  {(oi  - <o)t,  - s } 
n In 


n=l 


co 

V 

/ c (to  - to 

4 i n n 


to)  cos  {(to  - to)t  - e } 
n In 


n=l 


-E 


Cn(V  sin  {((V"  En^  cos  {(“-  u)t  - e } 
11  1 n n 1 n 


n=l 


0 


(A. A. 16) 


82 


E lxclxs2>  ‘ - E 


7 c (oj  - 
/ ■ n n 


u)  sin{(ui  - uj)t 
n 


1~  e } y c 
1 n / . n 


sin{(cu  - oj)  (t,+  x)-e 
n In 


n=l 


n=l 


= - y"1  4 c2  (a*  - u 
L i 2 n n 


n>)  COS  ( OJ  - 0))  T 
n 


n=l 


Thus,  E [x^x^l  equal  to  the  derivative  of  with  respect  to  t with 
opposite  sign.  Hence,  it  can  be  written  as 


E [xclxs2]  = 


(a)  - to)  S(w)  cos  (co  - o))t  do)  = - p 


1 


(A. A. 17) 


E [xclxs2'  ' 


- E 


v 

/ C (u)  - U) 

L I n n 


(jj)  sin  { ( oo  - 


w)t  - e } 
l n 


n=l 


uo 

/ c (u>  - 
/ ■ n n 


u)  cos  {(u)  - o>)  (t.  + x)  - e } 

n 1 n 


n=l 


* Z c» ( V 


-2 

w)  sin  ( oo  - o))t 
n 


n=l 


Thus,  E [xclxs2l  is  equal  to  the  derivative  of  -r^  with  respect  to  x. 
That  is, 


E ^xclxs2' 


/ 


-x  2 


(w  - u)  S ( oo)  sin  (oj  - co)  x dw  = 


(A. A. 18) 


83 


E [xc2xc21  ' - E 


t- 


cos  {(a)  - U))(t.+  t)  - s } 
n 1 n 


n=l 


oo 

EC  (ill  - 

n n 


(i>)  sin  {(u  - oj)  (t.+  t)  - e } 
n 1 n 


n=l 


-E 


c (a)  - oj)  COS  {(ill  - (D)(t.+  t)-  e } sin{  (oo  - uj)(t,+  x)-e  } 
n ii  n i n n 1 n 


n=l 


(A.A.19) 


E [xc2Xsl]  - E 


GO 

I, 


COS  { ( OJ  - oj)  (t 

n 


.+  t)  - e } ) C 

1 n Z_j  n 


sin  { (oo  - u)t  - e } 
n In 


n=l 


n=l 


oo 

--ifx 


sin  ( oo  - oj)t 
n 


n=l 


Thus,  E [XC2XS^]  equal  to  nQ  given  (A. A. 10)  with  opposite  sign.  That  is, 
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Thus,  E[xc2]  is  equal  to  the  second  moment  of  the  spectral  density 
function  about  its  mean  frequency,  denoted  by  y , That  is, 
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Thus , we  have 
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(A. A. 32) 
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Similarly  the  impressions  below  can  be  proved. 
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APPENDIX  B WAVE  SPECTRUM  COMPOSED  OF  TWO  SPECTRA  HAVING  SYMMETRIC  SHAPE 


In  order  to  simplify  the  covariance  matrix  given  in  Eq.  (3.5),  a 

given  wave  spectrum  is  decomposed  into  two  parts;  each  being  symmetric 

about  its  modal  frequency  as  illustrated  in  Figure  3.1.  One  part  contains 

primarily  the  lower  frequency  energy  components  of  the  wave  spectrum, 

while  the  other  contains  the  higher  frequencies  of  the  spectrum.  For 

this,  express  each  decomposed  spectrum  in  the  form  of  a normal  probability 

distribution  with  mean  w . and  variance  . That  is. 

mi  i ’ 


S.(m) 


/ 2tt  a. 

l 


exp  < 

L 


mi 


i = 1,  2 


(A.B.l) 


where,  mQi  is  the  area  under  the  spectrum  so  that  the  integration  of  Eq . 

(A.B.l)  over  the  frequency  domain  is  equal  to  m . . The  mean  value  to  is 

01  mi 

the  modal  frequency  of  S.(m).  By  letting  the  value  of  S.(o))  at  the  modal 
frequency  be  a the  unknown  variance  can  be  obtained  as 


(A.B.2) 


Thus,  the  symmetric  shape  can  be  expressed  by 


S.(u>) 


(A.B.3) 
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The  second  moment  about  the  mean  , of  the  spectrum  can  be  evalu- 
ated as 


y 


21 


(m  - 


. ) 

mi 


S.(m) 


dco 


(A.B.4) 


Two  conditions  are  considered  in  representing  the  shape  of  an  arbi- 
trarily given  spectrum  by  the  sum  of  two  spectra,  S^rn)  and  S 2 ( ou) , given 
in  Eq , (A.B.3).  These  are: 

(a)  The  area  under  the  original  spectrum,  mQ,  is  equal  to  the  sum 
of  mQj  and  mo2,  each  pertaining  to  the  area  under  the  spectrum 

S^(oj)  and  S2  (<*>),  respectively. 

(b)  Since  the  parameter  u2  defined  in  Eq . (3.6)  plays  a significant 
role  in  developing  the  probability  density  function  of  time 
duration  Ta+,  the  sum  of  y^  and  h22  evaluated  for  S^(rn)  and 
S2(u)  respectively,  is  maintained  constant  and  is  equal  to  y2 
of  the  original  spectrum. 

The  procedure  to  determine  the  parameters  of  the  two  decomposed  spec- 
tra is  as  follows: 

(1)  Divide  the  given  spectrum  S(oi)  into  two  parts  and  evaluate  the 
areas  mQl  and  mo2  (where  mol  + mo2  = mQ).  Choose  the  modal  frequencies 

“ml  and  “m2*  and  assume  al  as  the  value  of  the  spectrum  at  the  modal 
frequency  . All  these  values  are  appropriately  chosen  and  will  be  used 
as  the  initial  values  for  the  computations. 

(2)  Evaluate  y21  from  Eq . (A.B.4)  and  obtain  y 22  by  subtracting  y 
from  y^  which  is  calculated  for  the  given  spectrum. 

(3)  Evaluate  a2  from  Eq.  (A.B.4)  for  the  thusly  obtained  h22« 
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(4)  By  applying  a non-linear  least  squares  fitting  technique,  repeat 

the  above  computation  procedure  for  various  values  of  mQl,  aj,  Up  and 

such  that  the  difference  between  the  shape  of  the  original  spectrum 

2 

and  the  sum  of  two  symmetric  spectra,  expressed  by  | S(o>)  - E S (m)  | 2, 

i=l  1 

is  minimal. 


APPENDIX  C ELEMENTS  OF  COVARIANCE  MATRIX  FOR  WAVES  WITH 

SPECTRUM  OF  SYMMETRIC  SHAPE 


The  covariance  matrix  of  a random  vector  X associated  with  wave 
profile  is  given  in  Eq.  (3.5).  If  the  shape  of  a given  wave  spectrum  is 
represented  by  the  sum  of  two  symmetric  spectra,  then  the  covariance 
matrix  can  be  drastically  simplified  since  all  elements,  h , h^,  and  of 
the  covariance  matrix  become  zero.  The  proof  is  given  in  the  below: 

Since  the  shape  of  a given  wave  spectrum  is  represented  by  the  sum  of 
two  symmetric  spectra,  the  element  of  the  covariance  matrix  can  be  written 
as  the  sum  of  two  elements  of  the  covariance  matrices  evaluated  for  each 
decomposed  symmetric  spectrum.  For  example,  hQ=  nQl  + ho2,  where  n . ( i 
= 1 and  2)  is  given  by 


r 


noi= 


S^(to)  sin  (to 


to  . ) T doj 
mi 


(A.C.l) 


where  S^to)  - spectrum  with  symmetric  shape  given  in  Eq . (3.7) 


= a.e 


-IT 


oi 


2,  x2 

(„  - V) 


^ - modal  frequency  of  S.(to) 

= value  of  S^(io)  at  the  modal  frequency 
mol  = area  under  the  spectrum  S^to) 
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By  letting 


, i \ 2 

■fl  — = b.  and  w - a>  , = a 

moi ‘ 1 ml 


(A.C.2) 


we  have 


S.(a)  = a.e  bia 
l i 


(A.C.3) 


and  thereby 
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noi-  I s-|(a)  * sin  TCT  da  = aA  | e bia  sin  ta  • da  (A.C.4) 
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mi 


—a) 


mi 


Since  it  is  assumed  that  the  symmetric  spectrum  S.(a)  is  narrow- 

banded  and  its  energy  is  concentrated  in  the  neighborhood  of  its  mean 

frequency,  we  may  write  q . as 

oi 


00 

v 


( S.(a)  • sin  xa  da  = a.  J e biCT  • sin  t a • da 


00 


(A.C.5) 


The  integrand  is  the  product  of  a symmetric  and  an  unsymmetric  func- 
tions, hence  ^oi  e9ual  to  zero.  Thus,  it  can  be  proved  that 


n^  = ni+r|T  = 0 . 

o ol  o2 


In  a similar  fashion,  ^21  defined  as  follows  can  be  proved  to  be 
zero.  That  is, 
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(A.C.6) 
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Thus,  n2  = n21  + n22  becomes  zero. 
For  ti  we  may  write 


CO 


0 


• cos  (w  - to  . )t  ‘dm 
mi 


_oo 


a 


-b,  a 
e i 


cos 


xcr  da 


00  00 


[ 

2 

1 -b . oz 

r 2 

x / -b.o 

- e i • cos  to 

L 2bi 

+ 2b  1 e i sin  to  do 

L 

1 J J 

— CO  —00 


(A.C.7) 


The  first  term  of  Eq.  (A.C.7)  becomes  zero,  and  the  second  term  is 
also  equal  to  zero  as  shown  in  Eq . (A.C.5).  Thus,  and  thereby 

becomes  zero. 


APPENDIX  D 


DERIVATION  OF  JOINT 
f(  Xcl’  Xcr  Xc2’ 


PROBABILITY  DENSITY  FUNCTION 
Xc2’  Xsl’  Xsl’  Xs2’  Xs2) 


Derivation  of  the  joint  probability  density  function  of  the  random 
vector  X'  = (x^,  x^x^,  x^,  xgl>  xgl,  Xg2>  xg2)  given  Eq.  (3.9)  is  as 
follows.  The  joint  normal  probability  density  function  of  a n—dimensional 
random  vector  X can  be  written,  in  general,  as 


f(V  V •••*„>  ■ 7—772  -77172  - 1 - »>'  I'hx  - »»  (a.d.d 

(2ir ) I X I 

where  X’  = (x.  , x„,  x ) 

l z n 

V = vector  of  mean  values 
2 = covariance  matrix 

|E|  = determinant  of  the  covariance  matrix  E . 

For  the  present  problem,  the  covariance  matrix  of  the  random  vector 
given  in  Eq . (3.5)  is  drastically  simplified  to  that  given  in  Eq . (3.8)  by 
representing  the  shape  of  a given  wave  spectrum  by  the  sum  of  two  symmet- 
ric spectra.  We  may  write  each  element  of  the  covariance  matrix  of  the 
random  vector  E to  represent  the  sum  of  two  elements  given  in  Eq.  (3.8) 

evaluated  for  each  symmetric  spectrum.  Then,  the  covariance  matrix  E may 
be  written  as 


A 0 


E = 


0 A 


(A.D.2) 
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where 
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etc . 


(A.D.3) 


The  inverse  of  the  covariance  matrix  £ can  be  obtained  as 


(A.D.4) 


From  (A.D.2)  we  have  1=1  = |A|2-  Since  the  mean  values  of  the  random 
vector  are  zero  the  joint  probability  density  function  of  X can  be  written 
following  the  expression  given  in  Eq.  (A.D.l)  as 
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(A.D.5) 


where  |a|  = {(u2+  (y^)  - v2}  {(y^  v2)  (yQ  + vq)  - v2}  (A.D.6) 
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(A.D.8) 


The  elements  Mjj  in  the  above  matrix  can  be  evaluated  from  Eq, 
(A.D.3)  as  follows: 
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By  using  Eq.  (A.D.8),  K can  be  written  as 
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>2  >2  '2  '2 
<xcl  + xc2+  xsl  + xs2> 
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(A.D.10) 


APPENDIX  E DERIVATION  OF  JOINT  PROBABILITY  DENSITY  FUNCTION 

f(R1#  R2;  a) 


The  joint  probability  density  function  of  R^  R2>  0 for  a 

specified  level  a is  given  in  Eq.  (3.16).  That  is, 


2tt 


tiff 


f(R.,  R-,  e 8 ; a) *2 

2 (2*)3  |A| 


I 


-P., 
e d4> 


(A.E.l) 


where 
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(A.E.2) 


By  integrating  Eq . (A.E.l)  with  respect  to  O^and  02>  we  have 

f(R1»  r2»  a)  = 
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(A.E.3) 
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where 


Q = 
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2 M24^M22~  M24^  R1R2  C0S  ^ 

(A.E.4) 


and 


B1  ‘ - 2 a2  [M22  <M22Mir  lfU>  + <M22M13+  M14  ”24*  cos  * 

- <MuM24-  M22M^>  cos2  4 - M24(M13M24+  M24)  cos3  *J 

B2  - - 2 a [h22(H22«12-  M24Mu)  - M14(H22-  M24)  cos  ♦ 

- M24(M12M24“  M22M14}  cos2  *]  (A.E.5) 

In  order  to  evaluate  the  average  number  Na+(Ta+),  which  is  the 

numerator  of  Eq . (3.1),  the  joint  probability  density  function  f(R  ,R  ; a) 

given  in  Eq.  (A.E.3)  must  be  further  integrated  with  respect  to  and 
? 

R2*  Since  this  integration  is  extremely  difficult  for  the  form  presented 
in  Eq.  (A.E.3),  we  may  apply  the  following  technique  used  by  Rice  (1958) 
so  that  Eq.  (A.E.3)  will  be  reduced  in  a form  feasible  to  carry  out  the 

» t 

integration  with  respect  to  R^  and  R2<  That  is,  change  the  random  vari- 

11  i i 

ables  Rt  and  R2  to  Uj  and  U2  by  letting  R^  l^-  ^ and  R2  = U2~  k2  and 
choose  kx  and  k2  such  that  the  coefficients  of  U1  and  U2  in  the  resulting 
expression  of  the  term  Q in  Eq.  (A.E.3)  become  zero.  Thus,  we  have 
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B2+  2^M22_  M24^  ^k2M22+  klM24C0S  ^ 


= 0 


= 0 


(A.E.6) 


From  the  above  equation,  kj  and  k2  can  be  determined  as 

_ 11  [(M12M22-  M14M24>  + «12M24-  M14M22>  cos  ll 


kl  k2 


2 2 
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22  24 


(A.E.7) 


The  following  relationships  obtained  from  Eq.  (A.E.9)  are  also  used 
in  simplifying  Eqs.  (A.E.4)  and  (A.E.5): 
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(A.E.8) 


where  |A|  is  given  in  Eq.  (A.E.6).  That  is, 
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M12M22‘  M14M24  = {U2(V  0 " VQ}  VW  V0} 
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Let  us  write 
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where  |p|  < 1 following  the  definition  of  M22  and  M 


*24’ 


(A. 


Then,  from  Eqs.  (A.E.7)  through  (A.E.15),  Q given  in  Eq . (A.E.4) 
written  in  terms  of  Uj,  and  U2  as 
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(A.E.16) 
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where,  the  random  variables  (Uj,  U2)  and  ( R2>  have  the  following 

relationships : 
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u2=  r2+ 


k 


k = 


QLV)  ( 11  m c rh  — \)  ^ 


(A.E.17) 


From  Eq.  (A.E.17),  Q given  in  Eq.  (A.E.16)  can  be  expressed  in  terms 

i i 

of  and  R2,  and  thereby  the  joint  probability  density  function 

f t 

fCRj,  r2;  <*)  can  now  be  written  as 
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(A.E.18) 


° < R.  < 


< r2  < ° 


00 


and 


00 


APPENDIX  F SERIES  EXPANSION  OF  BIVARIATE  NORMAL  DISTRIBUTION 


The  joint  probability  density  function  of  two  normally  distributed 
random  variables  is  given,  in  general,  by 


f(x,y)  = 


2tt  0^2  (1-p2/2 


x exp 


2( 1-p2 ) 


(x  - v^)' 


- 2p 


(x  - y. ) (y  - 


t1 2 ) (y  - u2) 


2 v 


ai  a2 


(A.F.l) 


2 2 

where,  U1 , U2  are  the  mean  values  and  cr^ , o2  are  the  variances  of  x and  y, 
respectively.  p is  the  correlation  coefficient  between  x and  y.  For 
^2=  0,  the  joint  probability  density  function  becomes 


f(x,y)  = 


2tt  (1-p  )' 


W72 


exp 


2( 1-p2) 


iy ! - 2p  y y + d) 


- 


(A.F.2) 


Cramer  (1946)  shows  the  series  expansion  of  the  joint  probability 
density  function  given  in  Eq.  (A.F.2)  as  follows: 


(A.F.3) 
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Ill 


where,  $(  ) is  the  standardized  normal  cumulative  distribution  function. 

For  instance,  ) is  given  as 


x/ct, 


^ = 7^: 

1 / 2ir 


exp  { - 


} dt 


(A.F.4) 


By  taking  the  first  term  of  the  expansion  given  in  Eq.  (A.F.3)  and  by 
letting  o^=  a^=  1,  we  have 


_1 

2tt(1-P2  )i/2 

that  is, 


exp 


12  2 1 - x + y 

— — (x  - 2pxy  + y ) f ~ — e 2 


(_  2(l-p7 


1 r 2 2n 

(x  - 2pxy  + y ) 


2.  2 


2tt 


(A.F.5) 


2( 1-p  ) 


2 2 (A.F.6) 

2 1/,  - x + y 

(1-pV2  e 2 


This  is  the  approximate  formula  given  in  Eq.  (3.22). 


APPENDIX  G DERIVATION  OF  EXPECTED  NUMBER  N (x  ) 

a-  a-' 


The  expected  number  (per  unit  time)  of  the  envelope  downward  crossing 

' i 

with  velocity  R^  followed  by  the  upward  crossing  with  velocity  R2,  denoted 
by  N (x  ),  can  be  written  by 

u—  CC“ 


R2;  a)  dR^R2 


-“  < R1  < 0 

0 < r2  < » 


(A.G.l) 


The  procedure  to  obtain  Na_(x  ) is  the  same  as  that  for  the  deriva- 
tion of  the  expected  number  Na+(xa+)  shown  in  Eqs.  (3.2)  through  (3.29). 
However,  the  difference  takes  place  in  the  transformation  of  random  vari- 

t t 

ables  from  (R^ , R2)  to  (U,  V)  given  in  Eq.  (3.20).  That  is,  for  the 
derivation  of  Na_(xa_  ),  the  following  transformation  is  made: 


Then,  the  joint  probability  density  function  of  the  random  variables 
U and  V for  a specified  level  a becomes 
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2 TT 


x exp 


r - v cos  j>) 

oo 

2 2 

y - v 
o o 


1 


2 ( 1 — p ) 


(u2-  2p  u v + v2) 


By  applying  the  property  associated  with  the  bi-variate  normal  proba- 
bility distribution  given  in  Appendix  F to  the  2nd  term  of  the  exponential 
expression  in  Eq.  (A.G. 2),  the  joint  probability  density  function 
f(u,  v;  a)  approximately  becomes 


i i 

By  transformation  of  the  random  variables  R^,  and  R2  to  U and  V, 
respectively,  Eq . (A.G. 1)  can  be  expressed  in  terms  of  U and  V as  follows: 


M 

- L1oo 

N (x  ) = — 
a-  a-  2 2 

y - v 
o o 


00  00 

f f 

I | u v f(u,  v;  a 


) du  dv 


o o 


(A.G. 5) 
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Then,  from  Eqs . (A.G.3)  and  (A.G.4),  we  have 


(A.G.6) 


Since  the  terms  associated  with  U and  V are  symmetric  we  can  write 


CO  00 


| du  dv 


(A.G.7) 


By  letting 


= 2 y 


(A.G.8) 


the  integration  given  in  Eq.  (A.G.7)  can  be  obtained  as 


f -Cy  - 2 Y u) 

J ^ 


du  = 1 - /2ir 


Y e 


2Y 


V 


$ ( /2  y) 


0 
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where 


z 


(A.G.9) 


Hence,  Na_(ta_)  can  be  evaluated  by 


APPENDIX  H DERIVATION  OF  ELEMENTS  OF  COVARIANCE  MATRIX 
OF  THE  RANDOM  VECTOR  X'  = (xc,  xc,  xc,  xg , xg,  xg) 


The  covariance  matrix  of  the  random  vector  X' 
can  be  written  by 


= (xc,  xc>  xc,  xg,  xs 


E = 


E [x^] 
c 


E [x  x ] 
c c 


E [x  x ] E [x  ] 


c c 


E [x  x ] E [x  x ] 
sc  sc 


E [x  x ] 
c s 

t II 

E [ x x ] 
c s 


"2 

E[xZ] 

s 


(A.H.l) 


where  xc,  xg  are  expressed  by 


x (t)  = £ c cos[(u)  - m)t  - e 

c . n L n nJ 

n=i 


x (t)  = E c sinf(u)  - u)t  - el, 
s , n L n nJ 

n=l 


(A.H.2) 


en  are  t^ie  P^ases.  cn  are  the  amplitudes  which  are  related  with  the 
spectral  density  function  by 


I 2 

■7T  c =S(u))Au). 
In  n n 


(A.H.3) 


In  Appendix  A of  this  study,  the  elements  associated  with  displace- 


ment and  velocity  were  evaluated 


Those  are 
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and 


E[xc2' 

1 = do,  E[x£  xc 

E [x^  ] 

I 

= P0 , E [x  x 

c 

2 c s 

' 2 

E[xZ] 

s 

= y2  ‘ 

0,  E[x  x ] = 0,  E[x  x ] = 0, 
c s c s 

0,  E [x2  ] = Mq>  E[x  x ] = 0 (A.H.4) 

s s s 


The  derivation  of  elements  of  the  covariance  matrix  associated  with 
acceleration  is  given  in  this  Appendix  as  follows: 


E [x  x ] 
c c 


r 00 


-E 


E 


n=l 


c COSHO) 
n 1 n 


w)t  - 


00 

2 

e } E c (to  - to)  cos{(<o  - to) t 
n , n n 1 n 

n=l 


£ c2(m  - to)2  cos2{(<o  - to)t  - e 1 
n=l  n n n n' 


= - £ 
n=l 


1 2 

T c 
2 n 


£ 

n=l 


_ 2 

( 0)  - to)  S ( to  ) Ato 

n n n 


The  bar  denotes  the  average  with  respect  to  time  t.  By  letting  Aton 

ft 

small,  E[xc  x£]  is  equal  to  the  second  moment  of  the  spectrum  about  its 
mean  frequency  and  it  becomes 


E [x 


ii 

x 


CO 


(to  - 


to)2 


S ( to)  dto  = - u„ 


(A.H.5) 


If 


E [x  x ] 
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c cos  { (to  - <o)t  -e  } 
n n n 1 


E c (to  - to)  sm{(io  - co)t  -e  } 

, n n 1 n n JJ 

n=i 


(A.H.6 ) 
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t it 

E[x„  xl  = -E 
c c 


E c (to  - m) 
L , n n 
n=l 


sin{(w  - m)t  - e } 
L n nJ 


00 

2 

x E c ( co  - u)  cos  { ( to 
, n n L n 

n=l 
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1] 


= - E 
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2, 

C ( CO 

n n 
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{( 


to  - u>)t  - e 
n nJ 


sin{(a)  - m)t  - e ] 
1 n n ‘ 


= 0 


(A.H.7) 


_ 00 

» it  r 

E[x  x ] = E E c (to  - m)  sin{(m  - co)t  - e } 
c s L n n L n nJ 

n=l 


x I c (to  - a) ) sin{(oi  - w)t  — e }1 
. n n 1 n nJJ 

n=l 


E c (m  - w)  sin  { ( o>  - ui)t  - e } = E — c (m  - m) 
, n n 1 n nJ  , 2 n n 

n= 1 n= 1 


_ 3 

E (w  - u>)  S ( co  ) Aui  . 
, n n n 

n=l 


Thus,  E[xc  x ] is  equal  to  the  third  moment  of  spectrum  about  its 
mean  frequency  which  can  be  expressed  by 


1 II 

E [x  x ] 
c s 


( co  - m) 


3 


S(w) 


dm 


3 


"2  r 00  _ o _ 

E[x  ] = E I E c (oj  - u>)  cos{(u)  - m)t  - e } 

c u . n n L n nJ 

n=l 


(A.H.8) 
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x 


00 

£ c (ui  -to)  cos  { (to  - co)t  - e } 

, n n 1 n nJJ 

n=l 


= £ i c2(to  - co)4  = E (to  - to)4  S ( a)  ) A to 

, 2 n n . n n'  n 

n=l  n=l 


n 2 

Thus,  E[xc]  is  equal  to  the  fourth  moment  of  the  spectral  density 
function  about  its  mean  frequency.  That  is, 


"2 

E[x2] 


00 

■/ 


— 4 

(to  - to)  S(co)  dto  = . 


(A.H.9) 


E[xc  xg]  - E [ £ cn(un  - cos  {(to  - to)t  - e } 
n=l 


—v  2 
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00 

E cr.^u)n  - sin{((0  - to)t  — e } I = 0 (A.H.10) 

n “ n n J 
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»»  r __ 

E[x  x ] = -E  £ c sin{(to  - co)t  - e } 
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°°  2 1 
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n=l  n n n nJ 


£ sin2{(to  - co)t  -el 

n=l  n n n nJ 


r 1 2(  -N2 

£ c (to  - to) 

, 2 n n 
n=l 


2 

£ (to  - to)  S(to  ) Ato 
n=l  n n n 


Thus,  E [xsxg ] is  equal  to  the  negative  second  moment  of  the 
spectral  density  function  about  its  mean  frequency.  That  is, 
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E [x 
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v ■ -/ 


(m  - u>)  S(w)  dm  = -y. 


(A.H.ll) 
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I c (a)  - m)  cos{(m  - m)t  - e } 

. n n 1 n nJ 
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OO 

£ c (m  - m)2  sin{(m  - m)t  - e }]  = 0 (A.H.12) 

, n n 1 n n J 
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s L.nn  1 n nJ 

n=l 
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, n n L n n 
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n=l  n=l 


i»2 

Thus,  E[x^]  is  equal  to  the  fourth  moment  of  the  spectral  density 
function  about  its  mean  frequency  which  is  denoted  by 


"2 

EUJ] 


S(m)  dm  = y^ 


(A.H.13) 


APPENDIX  I DERIVATION  OF  JOINT  PROBABILITY  DENSITY  FUNCTION 
f(xc,  xc,  xc,  xg,  xs,  xs) 

The  derivation  of  the  joint  probability  density  function  of  the 

f II  f ?» 

random  vector  X'  = (xc,  xc,  xc,  xg,  xs,  xg)  given  in  Eq.  (4.7)  is  as 
follows : 

The  joint  normal  probability  density  function  of  a 6-dimensional 
random  vector  X with  zero  mean  can  be  written  as 


fCxj,  x2,  ...  x&)  = 


(2ir)' 


1/2 


exp 


{-T* 


,-l 


X} 


(A.I.l) 


where  X = (x^,  X2,  ...  x^) 

X = covariance  matrix 

|X|  = determinant  of  the  covariance  matrix 
For  the  present  problem,  the  covariance  matrix  X is  given  in 
Eq . (4.6),  and  we  can  write 


(A. I. 2) 


where 
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0 
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-u 
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-y 
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(A. I. 3) 
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The  inverse  of  the  covariance  matrix  £ can  be  obtained  as 


(A. I. 4) 


Since  the  determinants  of  the  two  submatrices  are  equal,  we  can 
write  |E|  = |a|^.  Hence,  the  joint  probability  density  function  of  X 
can  be  written  following  the  expression  given  in  Eq.  (A.I.l)  as 


f(xc, 


I If 

xs>  xc>  xs)  = 


(2ir)' 


|A| 


where, 
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s 
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= (xc,  xg , xc)  A 

X 
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X 

c 

tt 

X 

c 

II 

XS 

A-*  and 


are  given  as  follows: 
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M12 
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(A. I. 5) 


(A. I. 6) 


(A. I. 7) 
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11 

- m 

12 

M13\ 

12 

M 

22 

" M23  j 

13 

- M 

23 

M33/ 

(A. 1.8) 

The  elements  in  the  above  matrices  can  be  evaluated  from  Eq , 

(A. I. 3)  as  follows: 
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(A. I. 9) 
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23 


u u, 

o 3 


From  Eqs.  (A. I. 7)  through  (A. I. 9),  we  can  write  K as 


1 r 22  '2  '2  "2  "2 

K = [M. . (x  + x ) + M„_  (x  + x ) + M__  (x  + x ) 

jA|  L 11  c s 22  c s 33  c s 


T f 


If  f » II 


+ 2 M (x  x - x x ) + 2 M (x  x + x x ) + 2 M_„(x  x - x x )]. 
12  c s c s 13  c c s s 33  c s c s J 


(A. I. 10) 


APPENDIX  J DERIVATION  OF  JOINT  PROBABILITY  DENSITY  FUNCTION  f(£,8) 


I II  I 


The  joint  probability  density  function  of  R,  R,  R,  0 is  given  by 
Eq.  (4.13).  That  is, 


f(R,  R,  R,  0) 


(2tt)3/2(|a|M33)1/2 
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w 


(A.J.l) 


2 '2  2 ' 2 "2  '2' 2 '2' A 

where  P = MnR4  + M^R  + R^  0Z)  + M3;}(R  + 4 R 0 - R 0 


ii  * 2 2 ' " 2 * 2 
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II  ! I O I O I O 
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23' 


(A.J.2) 


By  multiplying  Eq . (A.J.l)  by  the  delta  function  <5(R)  and  by 

integrating  the  product  with  respect  to  R,  we  have 


+» 


f(R,  0,  R,  8) 


■/ 


6(R)  f(R,  R,  R,  9)  dR 
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(2tt)3/2(| A|  M33)1/2 
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2|Af 


(A.J.3) 


2 2 ' 2 "2  2 ' A it  12 

where  P'  = R + M22  R 0 + M3;j(R  + R 0 - 2 R R 0Z) 


22 


+ 2 M12  R2  0 + 2 M13(R  R - R2  02) 


124 


125 


tl  1 o t o 

+ 2 M23(R  r e - R QJ ) 

0 < R < °°,  < R < 0,  — 00  < 9 < °°  (A.J.4) 

By  using  the  joint  probability  density  function  given  in 
Eq.  (A.J.3),  we  can  evaluate  the  joint  probability  density  function 
given  in  Eq.  (4.19)  which  is  repeated  here;  namely, 


f(S,6) 


a c ae 


T {1 
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— 00 


Rf(5,  0, 


If 

R, 


f M 

0)  dR 


H If  f II  f 

Rf (R,  0,  R,  0)dR  dR  d0 


0 < 5 < °°,  ~°°  < 0 < 00 


For  evaluating  the  numerator  of  the  above  equation,  we  will 

II 

transform  the  random  variable  R to  H defined  as 


it 

R 


(lM) 1/2 

l m33J 


H 


(A.J.5) 


Then,  we  can  express  the  numerator  as 


Rf(5,  0, 


h f ( 5,  0, 


h, 


0)  dh 


(A.J.6) 


In  carrying  out  the  integration  involved  in  the  above  equation,  the 


following  formula  will  be  used: 
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00  2 
J h e"h2  + 2Cvh  dh  -i(l  + £v/i  eUv)  {1  + KSv)}).  (A.J.7) 


Thus,  the  numerator  of  Eq . (4.19)  becomes 
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/it  J 0 


(A.J.9) 


Next,  by  applying  the  integration  formula  given  in  Eq.  (A.J.7),  the 
denominator  of  Eq.  (4.19)  can  be  written  as 


4-00  OO  U 1/9 

r r C " « t « 1 |a| 

I I Rf(R,  0,  R,  a)  dR  dR  d9 L-1 372 
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(A.J.10) 
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Thus,  the  joint  probability  density  function  of  0 and  R can  be 
obtained  from  Eqs.  (A.J.8)  and  (A.J.10)  as 


(A.J.ll) 


APPENDIX  K PHASE  VELOCITY  k AND  TIME  INTERVAL  Ta 


In  this  Appendix,  the  phase  velocity  0 will  be  expressed 
approximately  in  terms  of  the  time  interval  between  two  successive  peaks 
of  the  envelope  above  a specified  level  a.  For  this,  the  wave  profile 
given  in  Eq.  (4.1)  may  be  written  as  follows: 


x(t)  = Re  E c e 

i n 

n=l 


i(w  t + e ) 
n n 


t,  i , N i0(t)  iwt 
= Re<p(t)  . e . e 


00  i(u>  - u>)t  + e 

■u  /\i^\t/  v n 

where,  P(t)  e = E c e 

n=l 


(A.K.l) 


Equation  (A.K.l)  implies  that  the  wave  profile  can  be  expressed  as 
a sinusoidal  carrier  with  the  amplitude  modulation  p(t)  ei^t^  varying 
slowly  with  time  in  comparison  to  the  carrier  frequency  (Longuet- 
Higgins,  1975). 

Let  us  consider  the  phase  difference  A0  between  two  successive 
envelope  peaks  both  of  which  are  above  a specified  level  a.  These  two 
peaks  are  separated  in  time  by  an  amount  approximately  equal  to  the  time 
interval  between  two  envelope  crossings  at  the  level  u,  denoted  by  Ta. 
Then,  we  may  write 


A0  = 0(t  + To)  - 0(t)  (A.K.2) 

By  expanding  the  right-hand  side  of  Eq.  (A.K.2)  into  a Taylor 
series,  we  have 
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A0  = 0 t 4 0 T2  + (A.K.3) 

a 2 a 

Since  0 ( t ) is  a slowly  varying  function  of  time,  the  higher  order 
derivatives  of  Eq.  (A. 4. 3)  may  be  neglected.  On  the  other  hand,  the 
phase  difference  between  two  successive  envelope  peaks  (not  necessarily 
above  the  level  a)  may  be  assumed  to  be  approximately  equal  to  2ir. 
Then,  the  phase  difference  between  two  successive  envelope  peaks,  both 
of  which  are  above  the  specified  level  a,  may  be  written  as  27rk.  Then, 
from  Eq.  (A.K.3),  the  following  relationship  can  be  derived: 

f 

0 T = 2itk  ( A.K.4) 

a 

The  unknown  parameter  k involved  in  the  above  equation  can  be 
determined  by  the  iteration  method  from  a consideration  of  the  proba- 
bility of  envelope  crossing  at  the  level  a.  Since  the  envelope  follows 
the  Rayleigh  probability  distribution,  the  probability  that  the  envelope 
R exceeds  a specified  level  a can  be  evaluated  by 


PrjEnvelope  exceeds  a} 


00 

■/ - 

J “0 


R 

2m0 


dR 


-mo 


= e 


(A.K.5) 


where,  mo  = S(m)  dm 


00 
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This  probability  is  equal  to  the  ratio  of  the  average  time  duration 
of  the  envelope  above  the  level  a,  denoted  by  ~a+,  and  the  average  time 
interval  between  two  successive  envelope  crossings,  x^.  Hence,  we  have 


e 


/x 


a 


(A.K.6) 


Ta+  can  be  evaluated  from  Eq.  (4.33)  of  Chapter  3 of  the  present 
study,  while  x^  can  be  obtained  from  Eq.  (4.24).  The  iteration  method 
is  applied  in  determining  the  k-value.  That  is,  the  probability  density 
function  f(xa|a)  given  in  Eq.  (4.24)  is  obtained  by  choosing  a k-value 
such  that  it  satisfies  the  relationship  given  in  Eq . (A.K.6). 
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